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Executive  Summary 


This  document  provides  a  description  of  Version  1.7  (revision  e)  of  the  Parameterized 
Peal-time  ionospheric  Specification  Model  (PRISM  1.7e),  a  global  data  assimilation  model 
intended  to  provide  near  real-time  ionospheric  specifications  for  space  weather  applications. 
PRISM  consists  of  a  semi-analytic  representation  of  a  parameterization  of  three  separate 
physically  based  computational  models  of  the  ionosphere.  PRISM  uses  both  ground  based  and 
space  based  data  to  modify  or  update  the  climatological  ion  density  profiles  to  produce  a  data 
driven  specification  of  the  state  of  the  ionosphere.  Currently,  PRISM  1.7e  is  operational  at  the 
Air  Force  Weather  Agency  (AFWA)  at  Offut  AFB  in  Nebraska,  and  is  also  a  component  of 
OpSEND,  an  operational  space  weather  software  suite  developed  by  the  Air  Force  Research 
Laboratory  (AFRL)  for  use  at  the  same  site.  In  operational  use  PRISM  typically  provides  either 
hourly  or  quarter  hourly  ionospheric  specifications  in  near  real  time.  PRISM  is  capable  of 
ingesting  and  using  electron  density  profile  (EDP)  parameters  (foFi,  hmF2,  foE,  and  hmE) 
produced  from  automated  true  height  analysis  of  digital  ionograms,  TEC  measurements  from 
GPS  receivers,  satellite  based  in  situ  plasma  measurements  (electron  density,  ion  composition, 
electron  and  ion  temperature,  and  ion  drift  velocity),  satellite  based  in  situ  observations  of 
precipitating  electrons  and  ions  in  the  high  latitude  regions,  as  well  as  ionospheric  parameters 
derived  from  observations  of  airglow  and  auroral  optical  emissions. 
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1.  INTRODUCTION 


We  have  modified  and  updated  the  Parameterized  Real-time  Ionospheric  Specification 
Model  (PRISM)  for  the  Air  Force  Research  Laboratory  (AFRL)  for  use  at  the  55th  Space  Weather 
Squadron  (55  SWXS)  and  by  the  Air  Force  Weather  Agency  (AFWA).  The  model  is  capable  of 
assimilating  both  ground  based  and  space  based  data  available  in  near  real  time  to  modify 
theoretical  climatology  provided  by  a  parameterization  of  several  regional,  first  principles 
ionospheric  models,  thus  providing  a  near  real-time  specification  of  the  ionosphere.  The 
theoretical  climatology  is  provided  by  a  composite  of  diumally  reproducible  runs  of  several 
physical  ionospheric  models:  (1)  the  Time  Dependent  Ionospheric  Model  (TDIM)  of  Utah  State 
University  (USU)  [Schunk,  1988],  (2)  the  low  latitude  F-region  model  (LOWLAT)  developed  by 
Anderson  [1973],  (3)  the  midlatitude  version  of  LOWLAT  (called  MIDLAT)  developed  by  D.  N. 
Anderson  and  modified  by  D.  T.  Decker,  and  (4)  an  E-region  local  chemistry  code  developed  by 
D.  T.  Decker  and  incorporating  photoelectrons  using  the  continuous  slowing  down  method 
[Jasperse,  1982]. 

PRISM  can  produce  either  regional  or  global  output.  The  output  grid  of  latitude  and 
longitude  (which  may  be  either  geographic  or  geomagnetic)  is  user  selectable.  (At  55  SWXS  the 
standard  PRISM  output  is  a  global  grid  with  1°  latitude  spacing  and  5°  longitude  spacing.) 
Profile  parameters  (f0F2,  hmF2,  TEC,  etc.),  or  complete  electron  density  profiles,  or  both  may  be 
output  to  the  specified  grid.  The  user  also  has  the  option  of  specifying  output  at  a  list  of  sites. 
Although  intended  for  validation  purposes  -  for  which  the  list  would  be  a  set  of  “ground  truth” 
sites  —  it  can  used  for  other  purposes,  such  as  approximating  the  ground  track  of  a  satellite. 

The  last  documentation  for  PRISM  described  Version  1.5.  The  current  version,  1.7e, 
corrects  a  number  of  bugs  that  were  detected  during  testing  at  AFRL  and  at  PRISM  has  been 
delivered  to  the  AFWA  and  is  presently  undergoing  transition  to  operational  code.  The  present 
documents  incorporates  some  of  the  material  in  the  previous  documentation  for  completeness, 
but  in  abbreviated  form,  and  concentrates  on  new  information  regarding  changes  to  the 
algorithms. 


1.1  Objectives 

The  primary  goal  of  PRISM  development  was  an  algorithm  for  assimilating  ionospheric 

data  in  near  real  time  in  order  to  produce  a  near  real  time  specification  of  the  global  ionosphere. 

The  data  to  be  assimilated  include 

(1)  bottomside  digital  soundings  from  the  Digital  Ionospheric  Sounding  System  (DISS), 

(2)  Total  Electron  Content  (TEC)  data  from  the  Ionospheric  Monitoring  System  (IMS)  and  other 
sources, 

(3)  in  situ  plasma  data  (densities,  temperatures,  and  drift  velocities)  from  the  SSIES  instrument 
on  DMSP  satellites, 

(4)  auroral  electron  and  ion  fluxes  from  the  SSJ/4  instrument  on  DMSP  satellites,  and 

(5)  electron  density  profile  information  deduced  from  observations  of  airglow  and  auroral  optical 
emissions  by  instruments  (SSUSI  and  SSULI)  expected  to  be  flown  on  future  DMSP 
satellites. 
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Because  the  SSUSI  and  SSULI  instruments  were  hardly  even  on  the  drawing  board  when  PRISM 
development  began,  no  explicit  input  specification  was  written  for  this  data.  Rather  it  was 
envisioned  that  these  instruments  would  provide  derived  electron  density  profile  (EDP) 
parameters  that  would  be  input  to  PRISM  as  if  they  came  from  DISS  or  SSIES  or  even  SSJ/4. 

The  need  for  a  global  specification  of  the  state  of  the  ionosphere  is  twofold.  First,  there 
are  operational  systems  that  need  to  correct  for  ionospheric  effects  in  real  time,  or  that  have 
operational  parameters  that  are  affected  by  the  ionosphere  and  must  be  adjusted  in  near  real  time. 
Second,  the  operation  of  many  systems  could  be  optimized  if  accurate  forecasts  of  ionospheric 
conditions  were  available  because  this  allows  the  operational  parameters  to  be  chosen  ahead  of 
time.  Any  ionospheric  forecast  algorithm  will  require  an  accurate  specification  of  the  current 
state  of  the  ionosphere  as  an  initial  condition,  which  PRISM  provides. 

In  addition  to  real  time  needs,  many  system  operators  need  post-event  analysis  to 
determine  whether  operational  problems  or  outages  were  caused  by  system  problems  or  by 
environmental  conditions.  PRISM  will  be  used  for  this  purpose  at  AFWA  as  well. 


1.2  Approach 

Ideally,  the  specification  of  the  current  state  of  the  ionosphere  would  be  obtained  directly 
from  real  time  observations  from  a  dense  network  of  satellite  and  ground  based  instruments. 
Unfortunately,  the  complexity  and  spatial  extent  of  the  ionosphere  precludes  the  deployment  of  a 
sufficiently  dense  network  of  observing  instruments.  Therefore,  any  practical  ionospheric 
specification  algorithm  must  be  based  on  an  ionospheric  model  with  parameters  that  can  be 
adjusted  on  the  basis  of  near  real  time  data.  Two  approaches  are  possible:  (1)  statistical  or 
empirical  climatological  models  or  (2)  numerical  simulations  based  on  physical  models,  (hi  this 
paper  “physical  model”  is  synonymous  with  “first  principles  numerical  model.”)  For  reasons 
described  below,  we  have  chosen  the  second  approach  (physical  models).  However,  practical 
considerations  (primarily  computational  speed)  dictate  that  the  algorithms  implemented  at  the 
AFWA  be  based  on  parameterized  versions  of  the  physical  models  (“theoretical  climatology”). 
Although  the  computational  requirements  are  less  stringent  than  when  PRISM  development  first 
began,  state-of-the-art  physical  models  still  require  high  speed  supercomputers  (especially  multi¬ 
processor  systems)  to  operate  faster  than  real  time. 

We  feel  strongly  that  a  comprehensive  physical  model  of  the  coupled  thermosphere  and 
ionosphere  can  produce  more  accurate  specifications  and  forecasts  than  can  statistical  or 
climatological  models.  The  causal  relationship  between  easily  monitored  solar  and  geophysical 
parameters  (e.g.,  Kp,  F10  7,  etc.)  and  a  particular  ionospheric  configuration  is  very  complex.  Any 
organization  of  historical  ionospheric  data  inevitably  averages  over  a  variety  of  ionospheric 
configurations  corresponding  to  similar  values  of  the  chosen  set  of  solar-geophysical  parameters 
(usually  only  one  or  two).  The  result  is  that  spatial  structure  tends  to  be  smeared  out  or  smoothed 
over,  and  the  resulting  model  is  unrepresentative  of  the  instantaneous  ionosphere.  If  a  physical 
model  contains  all  of  the  relevant  physics,  and  if  the  inputs  are  realistic,  then  it  will  produce  more 
realistic  representations  of  instantaneous  ionospheric  structure.  However,  there  is  a  difference 
between  a  realistic  representation  and  an  accurate  one. 
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In  order  to  accurately  simulate  a  time  dependent  phenomenon  like  the  thermosphere  and 
ionosphere,  a  physical  model  needs  an  accurate  specification  of  the  initial  conditions  and  an 
accurate  representation  of  the  energy  and  momentum  flux  at  the  boundaries.  For  the  purposes  of 
providing  a  specification  model,  it  is  the  energy  and  momentum  input  that  is  crucial.  If  the 
model  is  run  long  enough,  the  effects  of  the  initial  conditions  are  lost  and  the  present  state  of  the 
model  depends  only  on  the  recent  history  of  the  energy  and  momentum  input.  These  include  the 
solar  EUV  (the  primary  source  of  ionization  outside  the  auroral  zone),  high  latitude  heating  of 
the  thermosphere  (which  affects  the  global  circulation  of  the  thermosphere),  high  latitude 
convection,  and  low  latitude  dynamo  electric  fields.  While  the  temporal  and  spatial  resolution  of 
the  observations  of  these  quantities  are  expected  to  improve  in  the  future,  they  will  probably 
never  be  sufficient  to  allow  accurate  simulation  of  the  ionosphere  without  additional  data  in  the 
form  of  direct  measurements  of  ionospheric  parameters.  As  a  practical  matter,  ionospheric 
simulations  must  be,  and  will  remain,  iterative  in  nature:  The  energy  and  momentum  input 
parameters  are  adjusted  until  the  simulation  agrees  with  observations  of  ionospheric  parameters 
to  some  level  of  accuracy. 

There  is  a  superficial  similarity  between  our  approach  and  a  climatological  approach. 
The  difference,  however,  is  that  we  begin  with  a  more  realistic  representation  of  the  spatial 
structure  of  the  ionosphere  than  climatological  models  can  provide.  The  parameter  adjustment 
process  should  not  compromise  this  advantage.  In  the  future,  as  more  powerful  computers  and 
more  efficient  model  algorithms  become  available,  the  parameterized  models  can  be  replaced  by 
actual  physical  models  to  produce  more  accurate  specifications. 
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2.  THE  PHYSICAL  MODELS 


During  the  original  development  of  PRISM,  four  separate  physical  models  were  used  as 
die  basis  of  PRISM:  (1)  a  low  latitude  F  layer  model  (LOWLAT),  (2)  a  midlatitude  F  layer 
model  (MEDLAT),  (3)  a  combined  low  and  middle  latitude  E  layer  model  (ECSD),  and  (4)  a  high 
latitude  E  and  F  layer  model  (TDIM).  In  the  current  version  of  PRISM,  MIDLAT  has  been 
replaced  by  LOWLAT,  which  is  used  (separately)  for  both  the  low  latitude  and  midlatitude 
ionospheric  representations.  All  of  the  models  are  based  on  a  tilted  dipole  representation  of  the 
geomagnetic  field  and  a  corresponding  magnetic  coordinate  system.  (Hereafter,  “latitude”  means 
“magnetic  latitude”  unless  otherwise  noted.)  All  four  models  use  the  MSIS-86  neutral 
atmosphere  model  [Hedin,  1987].  Chemical  reaction  rates,  collision  frequencies,  and  similar 
data  are  consistent  among  all  the  models. 

2. 1  The  Low  and  Mid-Latitude  F  Layer  Model 

The  low  latitude  F  region  model  (LOWLAT)  was  originally  developed  by  Anderson, 
[1973].  (See  also  Moffett,  [1979]).  It  solves  the  diffusion  equation  for  0+  along  a  magnetic  flux 
tube.  Normally,  the  entire  flux  tube  is  calculated  with  chemical  equilibrium  boundary  conditions 
at  both  feet  of  the  flux  tube.  A  large  number  of  flux  tubes  must  be  calculated  in  order  to  build  up 
an  altitude  profile. 

Since  heat  transport  is  not  included  in  this  model,  ion  and  electron  temperature  models 
must  be  used.  For  the  PRISM  development  effort  we  chose  the  temperature  model  of  Brace  and 
Theis  [1981].  The  Horizontal  Wind  Model  (HWM)  of  Hedin  [1988]  was  used  to  describe 
thermospheric  winds.  We  used  production,  loss,  and  diffusion  rates  as  specified  by  Decker  et  al. 
[1994]. 

The  critical  feature  incorporated  in  the  low  latitude  model  is  the  dynamo  electric  field. 
The  horizontal  component  of  this  field  drives  upward  convection  through  ExB  drift,  and  this  can 
significantly  modify  profile  shapes  and  densities.  This  phenomenon  is  responsible  for  the 
equatorial  anomaly,  crests  in  ionization  on  either  side  of  the  magnetic  equator  at  ±15-20° 
magnetic  latitude.  In  the  current  version  of  PRISM  (Version  1.7e)  the  ExB  driven  vertical  drift 
used  for  these  calculations  was  based  on  the  empirical  models  derived  from  data  from  the 
Atmospheric  Explorer-E  (AE-E)  satellite  [Fejer  et  al.,  1995],  which  are  consistent  with  the  drifts 
measured  at  Jicamarca  [Fejer,  1981;  Fejer  et  al.,  1989]  but  include  longitudinal  variations  as 
well.  We  used  the  Fejer  et  al.  [1995]  empirical  drifts  for  moderate  and  high  solar  activity. 
Following  the  their  discussion,  we  modified  these  drifts  by  reducing  or  eliminating  the  pre¬ 
reversal  enhancement  for  low  solar  activity.  Horizontal  drifts  were  neglected  in  the  PRISM  runs. 
Of  course,  for  the  midlatitude  runs,  the  drift  is  irrelevant  (as  long  as  the  drift  vanishes  at  high 
altitudes  corresponding  the  flux  tubes  whose  feet  are  in  middle  latitudes). 

Since  its  original  development  this  model  has  undergone  extensive  validation  by 
comparison  with  data.  The  most  recent  such  comparison  is  Preble  et  al.  [1994]  using  electron 
density  profiles  measured  by  the  incoherent  scatter  radar  facility  at  Jicamarca,  Peru. 
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2.2  The  Low  and  Midlatitude  E  Layer  Model 

The  low  and  mid-  latitude  E  region  model  (ECSD)  was  developed  by  Dwight  T.  Decker 
and  John  R.  Jasperse  and  incorporates  photoelectrons  calculated  using  the  continuous  slowing 
down  (CSD)  approximation  [Jasperse,  1982].  Ion  concentrations  are  calculated  assuming  local 
chemical  equilibrium.  A  small  nighttime  source  is  included  to  ensure  that  an  E  layer  is 
maintained  throughout  the  night. 

2.3  The  High  Latitude  Model 

The  high  latitude  model  (incorporating  both  E  and  F  layers)  is  the  Utah  State  University 
(USU)  Time  Dependent  Ionospheric  Model  (TDIM).  (See  Schunk  [1988]  for  a  review.)  This 
model  is  similar  to  the  low  and  middle  latitude  models  except  that  the  flux  tubes  are  truncated 
and  a  flux  boundary  condition  is  applied  at  the  top.  In  addition,  the  flux  tubes  move  under  the 
influence  of  the  high  latitude  convection  electric  field.  In  the  low  latitudes,  because  the  magnetic 
field  is  mainly  horizontal,  the  effect  of  the  electric  field  is  primarily  to  move  the  ionization  in 
altitude.  In  contrast,  the  high  latitude  magnetic  field  is  mainly  vertical,  and  the  electric  field 
driven  convection  is  horizontal.  Like  LOWLAT,  this  model  has  a  long  history  and  has  been 
validated  by  numerous  comparisons  with  data. 

TDIM  includes  an  E-layer  model  that  incorporates  the  effects  of  ionization  by 
precipitating  auroral  particles.  The  ion  production  rates  used  were  calculated  using  the  B3C 
electron  transport  code  [Strickland,  1976;  Strickland  et  al.,  1994]  and  incident  electron  spectra 
representative  of  DMSP  SSJ/5  data.  The  characteristics  of  the  electron  spectra  were  taken  from 
the  Hardy  et  al  [1987]  electron  precipitation  model.  The  high  latitude  convection  patterns  were 
those  developed  by  Heppner  and  Maynard  [  1 987]  for  southward  directed  B. . 
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3.  PARAMETERIZATION  OF  THE  PHYSICAL  MODELS 


Parameterization  of  the  physical  models  proceeded  in  two  steps.  First,  the  models  were 
used  to  generate  a  number  of  output  "databases"  for  a  discrete  set  of  geophysical  conditions. 
Each  database  consists  of  ion  density  profiles  on  a  discrete  grid  of  latitudes  and  longitudes  for  a 
24  hour  period  in  UT.  Second,  to  reduce  storage  requirements,  the  databases  were  approximated 
with  semi-analytic  functions.  These  two  processes  are  described  in  the  following  subsections. 

3 . 1  Geophysical  Parameters 

All  the  physical  models  were  parameterized  in  terms  of  season  and  solar  activity.  The 
middle  and  high  latitude  models  were  also  parameterized  in  terms  of  magnetic  activity,  while  the 
high  latitude  model  was  additionally  parameterized  in  terms  of  the  sign  of  the  interplanetary 
magnetic  field  component  By.  (The  high  latitude  model  was  only  run  using  B.  southward. 

Northward  B.  conditions  are  modeled  using  the  low  magnetic  activity  databases.)  For  the  middle 
and  low  latitudes,  the  F  layer  (0+)  and  the  E  layer  (NO+and  Oj)  were  computed  and 
parameterized  separately.  The  high  latitude  model  (TDIM)  produced  all  three  ions  at  once. 

Due  to  time  and  computer  resource  limitations,  only  a  few  values  of  each  parameter  were 
used.  The  season  "values"  are  the  June  and  December  solstices  and  the  March  equinox  (which 
also  "stands  in"  for  the  September  equinox).  We  hope  to  change  from  seasonal  to  monthly 
values  in  a  future  version  of  PRISM.  The  values  of  the  other  parameters  are  summarized  for 
each  latitude  region  in  Table  1.  Note  that  the  USU  TDIM  and  LOWLAT  produce  output  in 
magnetic  local  time  (MLT),  while  MIDLAT  and  ECSD  produce  output  in  magnetic  longitude. 
Since  the  two  coordinates  are  readily  interconvertible,  we  will  ignore  the  distinction  and  refer 
only  to  magnetic  longitude  in  the  description  that  follows. 


Table  1 :  Geophysical  Parameter  Values 


Solar  Activity 
(^0.7) 

Magnetic 

Activity 
( Kp ) 

IMF  By 

direction 

Number  of 

databases 

Low  Latitude  F  layer 

70,  130,  210 

N/A 

3  6° 

Midlatitude  F  layer 

70, 130,210 

1 

N/A 

54* 

Low  &  Midlatitude  E  layer 

70, 130, 210 

EfSHM 

N/A 

54c 

High  Latitude  E  &  F  layer 

70, 130,210 

1,3.5,  6 

+.- 

324rf 

"3  seasons  x  3  solar  activities  x  4  longitude  sectors 

b3  seasons  x  3  solar  activities  x  3  magnetic  activities  x  2  hemispheres 

c3  seasons  x  3  solar  activities  x  3  magnetic  activities  x  2  species 

d3  seasons  x  3  solar  activities  x  3  magnetic  activities  x  2  By*s  x  3  species  x  2  hemispheres 
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3.2  Representation  of  the  Databases 

When  the  models  are  run  for  any  one  set  of  geophysical  parameters  (e.g.,  June,  FX01  = 
130,  Kp  =  1),  they  produce  ion  densities  (0+,  NO+,  and  07)  on  a  four  dimensional  grid. 
MIDLAT  and  ECSD  use  a  grid  of  magnetic  latitude  (A),  magnetic  longitude  ((p),  altitude(z),  and 
Universal  Time  (r ).  TDM  uses  magnetic  local  time  (MLT  or  y/)  instead  of  magnetic  longitude, 
while  LOWLAT  uses  MLT  instead  of  UT.  In  order  to  make  this  mass  of  numbers  more 
manageable,  we  produced  a  semi-analytical  representation  of  each  database.  The  space  and  time 
grid  parameters  are  summarized  for  each  latitude  region  in  Table  2. 

Due  to  the  computer  resource  requirements  of  the  low  latitude  F  layer  code,  it  was  used  to 
generate  databases  at  four  discrete  longitudes  (corresponding  to  longitude  sectors  for  which  ExB 
drift  measurements  were  available).  Each  longitude  sector  was  parameterized  separately,  and  the 
necessary  longitude  interpolation  is  carried  out  in  PM  and  PRISM  during  execution,  as 
described  below.) 

Because  we  were  trying  to  represent  discrete  data  (rather  than  continuous  functions),  and 
because  we  were  working  with  regional  rather  than  global  data  sets,  we  felt  that  the  usual 
spherical  harmonic  expansion  techniques  were  not  appropriate.  Instead  we  concentrated  on  the 
use  of  orthogonal  functions  of  discrete  variables. 


Table  2:  Horizontal  Grid  Parameters 


Latitude  Region 

Magnetic 

Latitude 

Magnetic 

Longitude 

UT 

Number  of 
altitude  profiles 
per  database 

Low  Latitude 

F  layer 

-44°  to  44°  in 

2°  steps 

30°,  149°,  250°, 
and  329° 

BH 

2,160 

Midlatitude 

F  layer 

34°  to  74°  and 
-34°  to  -74°  in 
4°  steps 

0°  to  345°  in 

15°  steps 

0100  to  2300  in 

2  hr  steps 

3,168 

Low  and 
Midlatitude 

E  layer 

-76°  to  76°  in 

4°  steps 

0°  to  345°  in 

15°  steps 

0100  to  2300  in 

2  hr  steps 

11,232 

High  Latitude  E 
&  F  layer 

51°  to  89°  and 
-51°  to  -89°  in 

2°  steps 

MLT:  0.5  to  23.5 
in  1  hour  steps 

0100  to  2300  in 

2  hr  steps 

5,760 

We  considered  the  use  of  modified  Chapman  functions  for  representing  altitude  profiles 
of  ion  densities.  These  functions  have  the  advantage  that  peak  height  and  peak  density  are 
explicit  parameters,  but  the  extremely  non-linear  nature  of  these  functions  necessitates  the  use  of 
non-linear  least  squares  fitting  methods.  While  such  methods  produced  excellent  representations 
of  individual  profiles,  the  variation  of  the  fitted  parameters  with  latitude,  longitude  (or  MLT), 
and  UT  was  unacceptably  noisy.  Consequently,  we  chose  to  use  Empirical  Orthomormal 
Functions  for  the  altitude  representation. 
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Empirical  Orthonormal  Functions  (EOFs)  have  been  used  extensively  to  represent 
meteorological  and  climatological  data  [Lorenz,  1956;  Kutzbach,  1967;  Davis,  1976;  and  Peixota 
and  Oort,  1991].  They  have  also  been  used  for  empirical  ionospheric  modeling  [Secan  and 
Tascione,  1984].  EOFs  are  described  in  Appendix  A.  They  have  the  advantage  of  providing  a 
representation  in  terms  of  linear  combinations  of  orthogonal  functions,  which  allows  for 
straightforward  determination  of  coefficients.  However,  because  peak  density  and  peak  height 
are  not  explicit  parameters  of  the  representation,  these  parameters  can  be  determined  only  by 
reconstructing  the  entire  profile  and  invoking  a  peak  finding  algorithm.  At  the  time  of  the  initial 
PRISM  development  effort  (1989-1991),  limitations  on  computer  recourses  and  other 
considerations  dictated  that  we  use  separate  sets  of  EOFs  for  each  ionospheric  database.  Since 
that  time  we  have  realized  that  we  could  have  accumulated  the  covariance  matrix  across  all  of  the 
databases  (as  long  as  we  used  the  same  altitude  grid  for  every  model  run)  and  thus  generate  a 
single  set  of  EOFs  for  the  entire  model.  We  have  since  done  this  (under  separate  funding)  for  a 
prototype  low  and  midlatitude  model  and  found  that  it  works  quite  well  with  only  a  small 
increase  in  the  number  of  EOFs  required  to  represent  the  ionospheric  profiles  to  high  accuracy. 
We  hope  to  be  able  to  reparameterize  the  physical  model  runs  using  a  single  set  of  EOFs  for  a 
future  version  of  PRISM.  This  would  permit  the  interpolation  on  EOF  coefficients  rather  than  on 
reconstructed  EDPs  as  is  required  by  the  current  version  of  PRISM. 

For  longitude  (or  local  time)  variations  (and  for  the  low  latitude  F  layer  UT  variation),  the 
obvious  choice  is  a  Fourier  series,  since  trigonometric  functions  retain  their  orthogonality 
properties  on  uniform  discrete  grids  and  because  the  data  is  periodic  in  the  independent  variable. 
These  worked  quite  well  for  the  high  latitude  models  under  all  conditions  and  for  the  low  and 
midlatitude  models  under  low  to  moderate  solar  activity  conditions.  However,  they  did  not  work 
well  for  the  low  and  midlatitude  models  under  high  solar  activity  conditions,  apparently  because 
the  EOF  coefficients  exhibited  exceptionally  large  gradients  at  dawn  and  dusk.  Therefore,  we 
decided  to  tabulate  the  coefficients  in  longitude  for  all  the  low  and  midlatitude  databases.  We 
have  experimented  with  using  EOFs  for  the  longitude  variations  in  our  prototype  low  and 
midlatitude  model,  and  may  use  this  approach  in  a  future  version  of  PRISM. 

For  the  latitude  variations,  we  chose  to  generate  grid-specific  orthogonal  polynomials 
using  the  algorithm  derived  in  Beckmann  [1973]  and  described  in  Appendix  B.  To  help  keep  the 
notation  straight,  we  summarize  it  in  Table  3. 


Table  3.  Notation  Summary 


variable 

index 

orthogonal  function 

index 

altitude 

h 

1  <  i  <  I 

EOF:  «„(*,) 

1  <m<M 

latitude 

\<j<J 

polynomial:  un{X) 

0<n<N 

longitude 

<Pk 

1  <k<K 

trigonometric0:  cos(/?<p)and  sin(jp^) 

0  <>p<P 

local  time 

1  <k<K 

trigonometric6:  cos(p^)and  sin(/?y/') 

0  <p<P 

UT 

1  <1<L 

trigonometric0:  cos(gr)and  sin(gr) 

0  <q<Q 

«Not  used  for  any  database  in  PRISM  1 .5,  but  may  be  used  in  future  versions. 
*Used  only  for  high  latitude  databases. 
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The  semianalytic  representation  of  each  database  was  generated  in  several  steps.  For  all 
ionospheric  regions,  the  first  step  was  the  determination  of  the  EOFs  from  the  ion  densities  in  the 
database  and  a  set  of  coefficients  for  representing  each  ion  density  profile  on  the 

latitude,  longitude,  UT  grid  (See  Appendix  A). 

M 

rroiM]  <la> 

171=1 

M 

[midlatitude  LOWLAT,  ECSD]  (lb) 

M 

[low  latitude  LOWLAT]  (lc) 

m=l 

where  z/5  X (pk,  y/k,  xp  and  y/t  are  all  points  on  the  model  output  grid,  and  g(^{zi)  is  the  m* 
EOF  evaluated  at  zr  (Note,  however,  that  a  different  set  of  (z(. )  functions  are  used  for  each 

ion,  for  each  set  of  geophysical  conditions,  and  for  each  model.) 

For  the  high  latitude  model  (TDIM,  both  E-  and  F-layers),  the  second  step  was  the 
generation  of  Fourier  coefficients  in  MLT,  a^(/ly,r;)  and  b{^(Xj,T,),  for  each  point  on  the 
latitude,  UT  grid. 

»,(Z„A;,^r,)*|;t{a«(A>,r1)coS(/.r)+i“(^,T,)Sm(p^)}g“(2:,)  [TOM]  (2) 

m= 1 p- 0 

For  the  low  and  midlatitude  models,  we  found  that  a  truncated  Fourier  series  often  introduced 
spurious  longitudinal  dependences,  apparently  driven  by  the  steep  gradients  at  dawn  and  dusk. 
The  effect  was  particularly  pronounced  at  high  solar  activity  when  the  day/night  contrast  is  the 
greatest.  Consequently,  for  these  models,  the  EOF  coefficients  remain  tabulated  in  longitude. 

For  all  models,  the  next  step  was  the  generation  of  orthogonal  polynomials  from  the 
latitude  grid  (Appendix  B).  For  the  high  latitude  model  (TDIM)  the  coefficients  are  and 

0-^np  (r , ) ,  and  the  ion  density  is  approximated  by 


MNP 


ns(Zi,X,y/,T,)«  / )  cos(p  y/)  +  P^np{r, )  sin(/?  y/)}gm  (z,>„  (X) 


m- 1  n- 0  p- 0 


[TDIM]  (3a) 


For  MTDT.AT  and  ECSD  the  coefficients  are  y{*l(Xk,  t,)  and  the  ion  density  is  approximated  by 
«,(zja,^,r/)*f;i^i(^,r/)gm(zI.K(A)  [midlatitude LOWLAT,  ECSD]  (3b) 

m= 1  n=0 

For  LOWLAT  the  coefficients  are  T](*l(<pk,  y/,)  and  the  ion  density  is  approximated  by 
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(3c) 


M  N 

ng(zi,A,Qk,ifrl)  tlow  latitude  LOWLAT] 

m= 1  n-0 

The  number  of  terms  in  each  series  are  listed  in  Table  4  for  each  region. 


Table  4:  Altitude  Grids  and  EPF's 


Database 

number  of 

altitude  points 

minimum 

altitude 

maximum 

altitude 

number  of 

EOF's 

low  latitude  0+ 

55 

160 

1600 

9 

midlatitude  0+ 

49 

125 

1600 

8 

low  &  midlatitude  NO+  &  02 

28 

90 

400 

7 

high  latitudeO+,  NO+,  &  Oj 

37 

100 

800 

6 

Note  that  in  none  of  these  cases  was  the  altitude  spacing  uniform. 


Because  of  the  extensive  use  of  tabulated  coefficients,  the  ion  density  at  an  arbitrary  point 
must  be  obtained  by  interpolation.  In  PIM  and  PRISM,  altitude  interpolation  is  quadratic,  while 
UT  interpolation  is  linear.  For  the  MIDLAT  databases,  the  longitude  interpolation  is  also  linear, 
as  is  the  local  time  interpolation  in  the  LOWLAT  databases.  However,  the  longitude 
interpolation  in  the  LOWLAT  databases  is  more  complicated.  First,  the  0+  profile  for  the 
desired  magnetic  latitude  and  local  time  is  reconstructed  for  each  of  the  four  longitude  sectors. 
Then  the  peak  height  and  peak  density  is  determined  for  each  profile.  The  peak  height  for  the 
desired  longitude  is  determined  by  Fourier  interpolation,  and  all  four  profiles  are  shifted  to  match 
the  interpolated  peak  height.  Then  Fourier  interpolation  is  used  at  each  altitude  to  obtain  the 
interpolated  ion  density  profile.  This  is  probably  not  ideal,  but  we  hope  to  replace  the  four 
longitude  sectors  with  12  or  even  24  in  a  future  version  of  PRISM. 

3.3  Merging  the  Regional  Models 

Because  we  used  four  different  regional  models  in  the  development  of  PRISM,  the 
models  must  be  merged  at  region  boundaries.  Specifically,  the  low  latitude  and  midlatitude  0+ 
models  have  to  be  merged  across  the  boundary  between  low  and  middle  latitudes,  while  all  three 
ions  (0+,  NO+,  and  0+2)  must  be  merged  across  the  boundary  between  midlatitudes  and  high 
latitudes. 

The  transition  from  low  latitude  0+  profiles  to  midlatitude  0+  profiles  takes  place 
between  34°  and  44°  in  both  hemispheres.  The  transition  is  accomplished  by  taking  a  weighted 
average  of  the  hmF2  values  from  the  two  models  in  which  the  weight  shifts  linearly  from  100% 
low  latitude  at  34°  to  100%  midlatitude  at  44°.  The  profiles  are  shifted  to  match  the  averaged 
hmF2  values  and  then  a  similar  weighted  average  of  the  shifted  profiles  is  taken  to  produce  the 
final  merged  profile.  Because  we  now  use  the  same  model  for  middle  and  low  latitudes,  the 
transition  is  smoother  than  in  earlier  versions.  No  transition  for  NO+  and  02  is  necessary  since  a 
single  model  was  used  for  these  ions. 
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The  transition  from  midlatitude  to  high  latitude  takes  place  over  an  8°  wide  zone  whose 
poleward  boundary  is  the  equatorward  boundary  of  the  trough.  The  transition  process  is  similar 
to  the  low  to  midlatitude  transition,  except  that  the  high  latitude  profiles  are  shifted  to  match  the 
hmF2  and  hmE  values  given  by  the  midlatitude  models.  The  final  profile  is  produced  by  a 
weighted  average  of  midlatitude  and  (shifted)  high  latitude  profiles. 

Although  PRISM  uses  geomagnetic  coordinates  internally,  it  can  produce  output  in  either 
geomagnetic  or  geographic  coordinates.  A  contour  map  of  NmF2  in  geographic  coordinates 
(cylindrical  projection)  for  1200  UT  on  21  June  2001  with  solar  activity  and  geomagnetic  activity 
indices  set  to  their  actual  values  is  displayed  in  Figure  1.  The  equatorial  anomaly  is  clearly 
visible  between  East  longitudes  30°  and  180°,  corresponding  to  local  times  of  1400  and  2400. 

PRISM:  NmF1:  2001.06.21,  1200  UT:  F1o.7=206.9,  SSNrff=  144.0,  Kp=  2.7 
-180  -135  -90  -45  0  45  90  135  180 


Figure  1.  Contours  of  NmF2  (in  units  of  105  cm"3)  in  cylindrical  equidistant  projection  from  PRISM  for 
1200  UT  on  21  June  2001.  The  equatorial  anomaly  is  clearly  evident  from  about  1400  to  2400 
local  time. 

The  high  latitude  region  is  more  clearly  seen  in  a  polar  projections  such  as  is  displayed  in 
Figures  2  and  3,  again  in  geographic  coordinates.  Figure  2  shows  the  Northern  (Summer) 
Hemisphere  and  the  Southern  (Winter)  Hemisphere  for  the  same  date  and  conditions  as  Figure  1. 
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PRISM:  N„FV  2001.06.21,  1200  UT:  F,a7-206.9,  SSN-f-1+4.0,  Kp-  2.7 


PRISM:  NmFa:  2001.06.21,  1200  UT:  F,o.7-206.9,  SSN„f-144A  Kp-  2.7 


4.  REAL  TIME  ADJUSTMENT  ALGORITHM 


The  Real  Time  Adjustment  (RTA)  algorithm  for  the  low  and  middle  latitude  regions  is 
different  from  the  algorithm  used  in  the  high  latitude  region.  This  is  partly  due  to  the  relatively 
complex  morphology  of  the  high  latitude  ionosphere  and  partly  due  to  an  evolution  in  our  ideas 
about  the  real  time  adjustment  process  during  the  development  of  PRISM. 

The  basic  philosophy  of  the  low  and  midlatitude  RTA  is  to  reproduce  the  input  data 
nearly  exactly  and  to  interpolate  correction  factors  smoothly  between  measurement  points.  Thus 
the  first  step  is  to  compare  PRISM  climatology  with  the  measurements  and  derive  corrections 
that  will  bring  PRISM  into  agreement  with  the  measurements  at  the  points  of  measurement. 
Then  the  next  step  is  to  calculate  an  “global”  correction  field  (or  an  interpolation  field)  that 
spreads  the  corrections  across  the  entire  low  and  midlatitude  region.  The  final  step  is  to  apply 
the  correction  field  to  PRISM  climatology  and  write  the  resulting  ionospheric  specification  to  the 
output  file.  (In  practice,  the  correction  field  is  calculated  “on  the  fly”  and  applied  to  each  output 
point  as  it  is  needed.) 

The  high  latitude  RTA  has  a  somewhat  different  philosophy.  Because  of  the  large 
amount  of  horizontal  structure  that  is  present  in  the  high  latitude  region  and  the  relative  paucity 
of  data  available  in  near  real  time,  the  RTA  involves  a  region-by-region  least  squares  approach. 
The  RTA  distinguishes  between  “indirect”  parameters  (such  as  an  effective  F10.7  or  Kp)  and 
“direct”  measurements  (such  as  a  measurement  of  electron  density  at  a  point).  It  treats  the 
indirect  parameters  as  adjustable,  and  varies  them  until  the  difference  between  PRISM  and  the 
direct  measurements  is  minimized  in  a  least  squares  sense.  In  practice,  the  RTA  uses  certain 
measurements  (e.g.,  energy  flux  of  particle  precipitation  and  horizontal  ion  drift)  to  set. 

It  should  be  noted  that  PRISM  was  originally  envisioned  as  assimilating  data  that  was 
nominally  obtained  at  the  instant  of  the  run  time.  That  is,  if  PRISM  is  run  for  1200  UT,  then 
only  data  taken  within  a  few  minutes  of  1200  UT  would  be  assimilated.  However,  it  quickly 
became  apparent  that  much  data  is  not  immediately  available  and  that  PRISM,  which  was  already 
anticipated  to  be  data  starved,  would  have  very  little  data  to  assimilate  if  operated  in  that  fashion. 
Therefore,  we  adopted  a  data  window  that  would  allow  PRISM  to  assimilate  all  data  taken  within 
a  four  hour  window  centered  on  the  nominal  run  time.  For  real  time  (or  near  real  time) 
operation,  this  effectively  means  that  PRISM  will  accept  data  that  is  up  to  two  hours  old,  while 
for  post  event  analysis,  PRISM  will  also  accept  data  that  was  taken  up  to  two  hours  after  the 
nominal  run  time. 

4.1  Available  Data 

The  near  real  time  data  available  for  use  in  the  adjustment  process  comes  from  the  Digital 
Ionospheric  Sounding  System  (DISS),  the  Ionospheric  Monitoring  System  (IMS),  and  a  suite  of 
Special  Sensors  on  the  DMSP  satellites.  The  DISS  network  consists  of  a  more  than  a  dozen 
digital  ionosondes  measuring  critical  frequencies,  critical  heights,  and  bottomside  profiles.  The 
DISS  ionosondes  are  derived  from  the  Digisonde  256  (D256)  developed  by  the  University  of 
Massachusetts  at  Lowell  (UML).  The  Ionospheric  Monitoring  System  consists  of  a  separate 
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network  of  dual  frequency  GPS  receivers  measuring  TEC.  In  addition  AFWA  receives  GPS 
TEC  data  in  near  real  time  from  a  network  of  approximately  24  GPS  receivers  that  are  part  of  the 
International  GPS  System  (IGS)  used  to  monitor  plate  tectonics  and  other  geophysical  data.  The 
DMSP  Special  Sensors  include  the  in  situ  plasma  properties  measured  by  the  SSIES  instrument 
and  the  precipitating  particle  measurements  of  the  SSJ/4  instrument.  In  the  near  future,  the  SSJ/4 
instrument  will  be  replaced  by  the  SSJ/5  instrument  and  a  new  set  of  ionospheric  remote  sensing 
instruments,  SSUSI  and  SSULI  will  be  added.  SSUSI  (a  multispectral  UV  imager)  and  SSULI  (a 
multispectral  limb  imager)  will  make  measurements  of  dayglow,  nightglow,  and  auroral  optical 
emissions  (mostly  ultraviolet).  The  observed  intensities  will  be  processed  on  the  ground  to 
deduce  ionospheric  properties  in  the  form  of  PRISM  profile  parameters  [e.g.,  Fox  et  al.,  1994]. 

Because  TEC  is  an  integral  quantity,  it  is  not  easy  to  incorporate  it  into  the  PRISM  real 
time  scheme.  In  PRISM  1.7,  we  have  chosen  to  convert  TEC  into  an  equivalent  point  datum. 
Our  method  for  doing  so  is  described  in  Appendix  C. 


4.2  Low  and  Midlatitude  Adjustment  Parameters 

The  PRISM  real  time  adjustment  algorithm  operates  on  six  parameters  that  prescribe  how 
an  electron  density  profile  is  to  be  modified  or  "corrected": 

1 .  tsf0F2 ,  the  correction  to  the  model  f0F2 , 

2.  Af0E,  the  correction  to  the  model  faE, 

3.  A hmF2 ,  the  correction  to  the  model  hmF2 , 

4.  A hmE,  the  correction  to  the  model  hmE, 

5.  ANwp,  the  correction  to  the  0+  density  at  a  specific  altitude  (i.e.,  the  DMSP  altitude) 

6.  A Htop,  the  correction  to  the  0+  scale  height  at  a  specific  altitude  (i.e.,  the  DMSP 
altitude) 


The  nominal  value  for  each  of  these  parameters  is  zero.  A  positive  (negative)  value  means  that 
the  model  value  must  be  increased  (decreased).  Using  the  available  near  real  time  data,  the  real 
time  adjustment  process  will  assign  non-zero  values  at  each  location  where  data  is  available  (the 
driver  sites). 

Parameters  1-4  are  based  on  direct  measurements  by  DISS  digital  ionosondes.  Parameter 
5  is  based  on  the  direct  measurement  of  the  0+  density  by  the  SSIES  instrument  on  board  the 
DMSP  satellite.  In  contrast,  Parameter  6  must  be  inferred  from  the  electron  and  ion  temperatures 
(Te,  Tt)  measured  by  the  SSIES  instrument.  At  midlatitudes,  PRISM  assumes  that  the  topside 
ionosphere  is  in  diffusive  equilibrium  and  calculates  the  topside  scale  height  from 

(4) 

The  corresponding  model  value  is  obtained  from  the  model  0+  densities  (n2,  n,)  at  the  grid 
altitudes  (z2,  zj)  immediately  above  and  below  the  DMSP  altitude: 

H  (5) 

]n(njn2) 

At  low  latitudes,  the  topside  scale  height  is  determined  by  diffusive  processes  alone,  so  the  scale 
height  cannot  be  inferred  from  the  SSIES  temperature  data. 


densities  (n2,  nx)  at  the  grid 
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After  PRISM  reads  all  of  the  input  data  files  and  writes  them  into  a  single  scratch  file,  the 
data  ingestion  process  proceeds  with  five  passes  through  the  combined  input  data  file.  The  first 
pass  searches  for  layer  height  data  (Parameters  3  and  4)  and  calculates  the  height  correction  for 
each  such  datum.  On  the  second  pass,  PRISM  searches  for  peak  density  data  (Parameters  1  and 
2)  applying  the  height  corrections  as  it  goes  along.  On  the  third  pass,  the  code  searches  for  TEC 
data  that  is  farther  than  the  parameter  DCL  (decorrelation  length,  currently  set  to  1000  km  at 
compile  time)  from  the  nearest  ionosonde  station  and  calculates  peak  density  corrections  as 
described  in  Section  4.5  below.  On  the  fourth  pass,  the  code  searches  for  SSIES  in  situ  plasma 
data  and  calculates  topside  corrections  (Parameters  5  and  6).  On  the  fifth  and  final  pass,  the  code 
searches  for  TEC  data  that  is  closer  than  DCL  to  the  nearest  ionosonde  but  farther  than  DCL 
from  the  nearest  SSIES  datum.  Any  TEC  data  that  meet  this  criteria  are  used  to  calculate  topside 
corrections  as  described  in  Section  4.5  below.  (Any  TEC  data  not  meeting  the  criteria  of  the 
third  or  fifth  pass  are  ignored.) 


4.3  Real  Time  Adjustment  of  the  Low  and  Midlatitude  Profile  Parameters 

In  PRISM  1 .7,  after  the  values  of  the  profile  adjustment  parameters  have  been  determined 
at  each  driver  site,  the  global  correction  field  is  determined  using  a  weighted  average  method. 

1.  Given  N driver  locations  and  the  associated  geomagnetic  coordinates  [X „,0„),  let  the 
unadjusted  PRISM  value  at  the  nth  point  be  un  and  the  measured  value  be  v„ . 

2.  The  correction  to  be  applied  to  the  unadjusted  PRISM  (i.e.,  PIM)  value  at  the  nth 
point  is  cn=vn-un. 

3.  At  any  other  point,  [X,  </>) ,  the  correction  to  be  applied  to  the  PIM  value  is 


EW»(A>^K 

c(A,^)  =  -=fc -  (6) 


where  the  wn(X,<jj)  are  weight  functions  that  depend  on  the  distance  measure  dn(X,j>) 
between  the  point  {X,<j>)  and  the  n,h  point(/l n,<f>n). 


dJIJ) 


N 


n'.M 


(7) 


In  PRISM  1 .7  we  have  used  the  following  distance  measure. 

(8) 

where  yn(X,<fi)  is  the  great  circle  distance  between  {X,<£)  and  {Xn,(f>n)'. 
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cos  yn  (A ,  <f>)  =  sin  A  sin  An  +  cos  A  cos  An  cos(^  -  <pn ) 


(9) 


As  a  practical  matter,  the  actual  weighting  function  used  in  PRISM  is  wn{A,(p)-\l  dn{A,(p) 
unless  dm(A,<p)  <  S  for  some  driver  station  m.  In  that  case,  wm  =  1  and  wn  =  0  for  n  *  m.  In 
PRISM  1.7,  S  is  set  to  1(T6.  It  should  be  noted  that  in  the  case  where  a  station  has  supplied 
several  measurements  during  PRISM’s  four  hour  data  window,  the  above  algorithm  effectively 
averages  the  corrections  for  that  station.  Consequently,  the  RTA  may  not  precisely  reproduce  a 
measurement  taken  at  the  nominal  time  of  the  run. 

Because  the  decorrelation  length  of  the  ionosphere  is  of  the  order  of  1000  km  or  less,  no 
interpolation  scheme  can  hope  to  accurately  reproduce  ionospheric  parameters  where  there  is  no 
data.  Clearly,  the  denser  the  data  net,  the  better  the  model  will  do.  Unfortunately,  the  data  sets 
available  so  far  are  rather  sparse,  particularly  over  the  oceans.  PRISM  provides  the  option  of 
introducing  “phantom  stations”  in  regions  with  no  data  in  order  to  force  PRISM  to  relax  to 
climatology  in  those  regions.  We  recommend  that  under  normal  circumstances  PRISM  be  run 
without  PHANTOM  stations.  (The  option  is  a  compile  time  option,  requiring  that  PRISM  —  or  at 
least  the  main  program  -  be  recompiled  in  order  to  exercise  this  option.)  For  completeness,  we 
describe  the  procedure  the  PRISM  follows  for  placing  phantom  stations. 

The  low  and  midlatitude  region  is  divided  into  32  rectangular  subregions  with  boundaries 
defined  by  the  following  parallels  of  latitude  and  meridians  of  longitude: 

parallels  of  latitude:  60°S,  30°S,  0°,  30°N,  and  60°N 

meridians  of  longitude:  30°E,  75°E,  120°E,  165°E,  210°E,  255°E,  300°E,  and 

345°E 

There  is  a  phantom  station  located  at  the  center  of  each  rectangle.  During  the  data  ingestion 
process,  PRISM  keeps  track  of  how  many  stations  are  located  in  each  rectangular  subregion,  and 
how  many  have  f0F2  data.  As  long  as  at  least  one  station  in  a  subregion  reports  a  value  for  f0F2, 
then  the  phantom  station  for  that  subregion  is  ignored.  If,  however,  no  station  in  the  subregion 
reports  a  value  for  f0F2,  then  the  phantom  station  is  assigned  a  value  provided  by  PRISM 
climatology  at  that  point.  A  similar  procedure  is  followed  for  hmF2,  F0E ,  and  hmE  data.  These 
phantom  stations  force  the  correction  field  to  relax  to  values  near  zero  in  regions  where  there  is 
no  data. 

4.4  Modifying  the  Low  and  Midlatitude  Model  Profiles 

The  interpolation  method  described  above  provides  a  global  correction  field  that  can  be 
used  to  calculate  the  profile  correction  parameters  at  any  location.  In  this  section,  we  describe 
the  way  in  which  the  profiles  are  adjusted  using  the  profile  correction  parameters.  First  the  layer 
heights  are  adjusted,  then  the  layer  peak  densities,  and  finally  the  topside  correction  is  applied. 

The  layer  height  correction,  using  parameters  3  and  4  is  a  simple  shifting  of  the  profiles. 
If  the  F  layer  correction  is  A hmF2 ,  then  the  0+  profile  is  shifted  so  that 

»„.”'(z)  =  V““(z-AA.Fa)  0°) 

A  similar  shift  is  applied  to  the  molecular  ion  profiles,  except  that  the  altitude  shift  is  A hmE. 
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Once  the  altitude  corrections  have  been  applied,  parameters  1  and  2,  the  critical 
frequencies,  will  be  used  to  scale  the  ion  density  profiles.  First,  the  peak  densities  of  the  E-  and 
F-layers  are  first  converted  to  critical  frequencies.  The  additive  correction  is  applied  to  each 
critical  frequency,  which  is  then  converted  back  to  a  density.  The  ratios  of  the  critical  corrected 
to  climatological  peak  densities  are  use  to  scale  the  molecular  and  atomic  ion  density  profiles. 


Rf 


NmFl 


(corrected ) 


KF2 


(climo) 


jy  £( corrected ) 

r>  _  T  tn 

E  t  *-,(climo) 


(11) 


The  ratio,  Rf,  is  used  to  scale  the  ion  density  profiles  (molecular  and  atomic)  above  the  F- 
layer  peak.  The  ratio,  Re,  is  used  to  scale  ion  density  profiles  (molecular  and  atomic)  below  the 
F-layer  peak.  Between  the  two  layer  peaks,  the  scaling  process  is  more  complicated.  First  we 
define  a  transition  altitude  htF 2  which  is  the  larger  of  the  (a)  altitude  at  which  the  0+  density  has 
fallen  to  e2N„F2  and  (b)  =  hmF2  -f  (hmF2-hmE)  =  hmE+\(hmF2  - hmE ) .  Then  we  define  a 


scale  factor  that  smoothly  varies  between  RE  at  hmE  and  RF  at  hiFi-  Explicitly,  the  density  scale 
factor  (applied  to  both  molecular  and  atomic  ions)  is 


z  >  hbF2 


S(z)  = 


EE  +  2  (  EF  Ee  ) 
Re, 


1-cos 


g) 

KKF2~KE  yj 


hmE  <z<  hbF2 
z,<  hE 

tn 


(12) 


At  the  altitude  boundaries  hmE  and  htFi,  S(z )  has  the  desired  properties 


S{h.E)  =  Rc 
dS\ 


dz 


=  0 


I  hmE 


and 


S(hbF,)  =  RF 


so  that  it  varies  smoothly  across  the  altitude  boundaries. 


(13a) 


(13b) 


Parameters  5  and  6  are  used  to  correct  the  topside  profile  based  on  ne,  Te,  and  Tx 
measurements  from  SSIES  on  DMSP  (nominally  840  km).  Let  NJz)  be  the  model  0+  profile 
( after  fJF2  and  hJF2  corrections  have  been  applied).  Further,  let  zp  =hmF2,  and  let  ztop=  the 
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altitude  at  which  A Ntop  and  A Htop  were  measured  (usually  the  DMSP  altitude).  Let  N£z)  be  the 
corrected  profile  based  on  ANtop  and  A Htop.  Nc(z)  must  satisfy  the  following  conditions: 


and 


Nc(z!op)  =  Ntop  =  Nm(ztop)+ANlop 


'  1  dNc 
\Nc  & 


'‘'top 


_1 _  1 

^top  Hm  (zlop  ^  +  A Htop 


(14a) 

(14b) 


In  PRISM  the  corrected  profile  is  obtained  from  the  model  profile  by  scaling  the  altitude  so  that 


where 


where 


and 


Nc(z)=Nm(() 


<•(*)- 


z, 


z<z„ 


(z-z  )2 

Zp  +  [2*1  -  ]{Z -  Zp  )  ■ +  K  -  *l]  Z  <  ztcp 


Z  —2 

top 


zx+RH{z-z[op), 


z>z, 


top 


Nm{zx)^Ntop 

z\~zp 

*i= - — 

Ztop  Zp 

.  _#.(*,) 

R H - S - 

**  top 

(  1  CN„\ 


Hm{z\)  V K  &  J2=zi 


(15a) 


(15b) 

(16a) 

(16b) 

(16c) 

(16d) 


Direct  substitution  of  z-ztop  in  Equation  (15)  verifies  that  £(ztop)  =  Zj  so  that  Equation  (14a)  is 
satisfied.  That  Equation  (14b)  is  satisfied  may  be  seen  from 


h.K) 


'  1  SN,' 

N  dz 

\  Jztop 

(  1 

Nm(t)  ar  <?z 


Jztop 


(17) 


H, 


top 


This  form  ensures  that  the  correction  and  its  first  derivative  vanish  at  the  peak  and  that  the  0+ 
concentration  and  its  slope  are  as  specified  at  ztop. 
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When  there  is  no  topside  scale  height  data  (as  at  low  latitudes),  PRISM  sets  Htop  to  the 

value 

Hwp  =  —Hm(zx)  (when  Htop  is  not  determined  from  data)  (18) 

x, 

so  that  the  altitude  scaling  becomes  linear 

£(z)  =  zp  +  x1(z-zp)  (when  Htop  is  not  determined  from  data)  (19) 

Now  let  us  examine  the  behavior  of  C,  (z) .  When  Htop  is  not  determined  from  the  data, 
£  is  clearly  well  behaved:  £  (z2 )  >  £  (z, )  if  z1  >  z,  or  in  other  words, 

^>0  (20) 

dz 


as  long  as  zx>  zp,  which  it  should  always  be.  On  the  other  hand,  in  the  quadratic  case, 

1,  2<zp 

(2x, -Rh)  +  2(Rh -xt)  Z  Zp  ,  zp<z<ztop  (21) 

2 top  Zp 

R-H>  2  ~  2 lop 

We  require  d<£ / dz  >0  for  zp<z<  ztop .  This  means  that  za ,  the  altitude  at  which  d£  /dz  =  0 , 
must  either  be  below  zp  or  above  zlop.  If  it  is  below  zp ,  then  we  require  that  it  represent  a 
minimum,  i.e.  that  d2£  /dz2  >0.  On  the  other  hand,  if  it  is  above  ztop ,  then  we  require  that  it 
represent  a  maximum,  i.e.,  that  d2C,  jdz2  <0. 


Ze*=Zp+\{Z,op-Zp) 

2x,  -  Rh 
x,  —RH 

(22) 

d __  2{^h  ~~x\ ) 
dz2  z  —z 

*top  p 

z>zp 

(23) 

Since  ztop  >  zp  is  almost  always  true,  z^  <  zp  only  if 

2x1  —  Rh  ^  q 

Xj  —Rff 

(24) 
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which  implies  that  \RH  <xx<R„,  since  RH>  0  by  definition.  On  the  other  hand,  >  ztop 
implies  that 


or 


2X)  (25) 

X\  —RH 

—5—  >1  (26) 

xi  ~Rh 


which  is  always  true  as  long  as  x,  >  RH  >  0 . 

As  an  example,  let  us  consider  the  following  case: 

Nm(z)  =  Nm(zp)ex p 


1-^ 

Hn 


-exp 


f  z—z  ^ 

p 

’  Hn 


(27) 


where  zp  and  Nm  have  already  been  adjusted  to  fit  the  data.  Suppose  that  we  have  measurements 
of  Ntop  and  Htop  at  ztop  =  zp  +  8 H0 .  Then  z\  is  determined  from 


Ntop  =  Nm  (Zl)~  Nm  (Zp)eXP)  2 


!_£lZ5l 

Hn 


-exp 


f  Z\-ZP^ 


H, 


o  )J 


or 


Zj  =  zp+H0 


l-21n 


N, 


top 


Mz>) 

Solving  iteratively,  we  have  the  “zeroth”  approximation 

z}°>  =*,+#„• 


-exp 


Hn 


AT 

l-21n 

top 

/  \ 

> 

_N-{zr)_ 

and  the  first  approximation 


l-21n 


N. 


top 


-exp 


-l  +  21n 


f  N  N 

v  top 


or 


Z?  =Zp+H0 


l-21n 


TV, 


top 


«.(*>) 


N, 


\2 


J-i 


top 


Mzp) 


(28) 


(29) 


(30) 


(31a) 


(31b) 
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and 


.  z'-z» ..  H,  , 

l-21n 

1 

•§ 

1 

-p~1 

f  Ntop  ] 

-2 " 

^top  ~  Z p  Z(op  —  Zp 

_Mzp)_ 

—  (' 

Nm(zn) 

'  p) ) 

(32) 


Suppose  Ntop  =  e~2Nm  (zp ) ,  then  z,(l)  =  zp  +  H0  |l  +  4  -  <?~ler4 } «  zp  +  4.99H0  and  x, 
With  zlop  =zp+  SH0 ,  x,  *  0.624 . 

Now 

1  1  dN 


so 


H.M  KM  & 
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(33) 


(34) 


£(z)  = 
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JU 

■8 
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V  P ) 
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-- 0.624 


top 


Z-Z„ 
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(36) 


Z>Z, 


/op 


An  example  of  the  midlatitude  profile  correction  algorithm  applied  to  a  single  profile  is 
shown  in  Figure  3,  for  which  we  have  used  Incoherent  Scatter  Radar  (ISR)  data  from  Arecibo  to 
simulate  simultaneous  Digisonde  and  SSIES  data.  In  this  case,  the  data  is  from  1800  UT  on  4 
October  1989.  In  Figure  3a,  the  ISR  data  is  compared  with  PRISM  climatology.  A  profile  from 
the  Ionospheric  Conductivity  and  Electron  Density  (ICED)  model  [Tascione  et  ai,  1988]  is 
shown  for  comparison.  The  actual  NmF2  is  a  factor  of  two  higher,  and  the  actual  hmF2  is  about  40 
km  higher,  than  the  climatological  values.  In  Figure  3b,  the  PRISM  profile  has  been  adjusted 
only  to  match  the  actual  NmF2.  In  Figure  3c,  both  NmF2  and  hmF2  have  been  adjusted  to  match 

the  actual  data.  Finally,  in  Figure  3d,  both  the  peak  and  topside  adjustments  have  been  made.  If 
we  define  an  RMS  density  error  as 
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then  PRISM's  RMS  density  error 
declines  from  46%  with  no  profile 
adjustment  to  25%  with  an  NmF2 
adjustment  only,  to  18%  with  both 
NmF2  and  hmF2  adjusted,  and  to  16% 
with  both  peak  and  topside  adjustments. 

Note  that  in  operational  use 
there  are  normally  many  data  sites,  and 
the  topside  and  bottomside  data  are 
seldom  co-located.  Operationally, 
PRISM  uses  all  available  hmF2  and  hmE 
data  to  establish  global  layer  height 
correction  fields,  then  it  uses  the 
available  NmF2  and  NmE  data  to 
calculate  global  profile  scaling  fields, 
and  finally  it  uses  the  topside  data  to 
calculate  global  topside  correction 
fields.  These  correction  fields  are  then 
applied  to  each  profile  on  the  output 
grid  as  described  above. 
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Figure  3.  The  PRISM  profile  adjustment  procedure 
illustrated  using  Incoherent  Scatter  Radar  (ISR) 
data  from  Arecibo.  (a)  No  adjustment:  The  RMS 
density  error  is  46%.  (b)  NmF2  adjustment  only: 

The  RMS  density  error  is  25%.  (c)  NmF2  and 
hmF2  adjustment:  The  RMS  density  error  is  18%. 
(d)  NmF2,  hmF2,  and  topside  adjustments:  The 
RMS  density  error  is  16%. 


4.5  The  use  of  TEC  data  in  the  low  and  midlatitude  RTA 

In  order  to  use  TEC  data,  PRISM  must  convert  the  TEC  measurement  into  an  equivalent 
point  measurement.  TEC  data  is  ingested  as  vertical  equivalent  TEC  at  the  Ionospheric 
Intersection  Point  (HP).  If  the  IIP  is  farther  than  DCL  (currently  1000  km)  from  the  nearest  DISS 
site,  then  the  TEC  is  converted  into  equivalent  values  of  NmF2  and  NmE  by  scaling  the  entire 

electron  density  profile: 
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The  corrected  peak  densities  are  converted  into  critical  frequency  parameters  and  critical 
frequency  corrections  are  calculated  as  if  the  data  were  obtained  from  ionosondes: 


1] 

V.B  = 


(39) 


From  this  point,  PRISM  treats  TEC  data  as  if  it  were  ionosonde  data  with  hmF2  missing. 

If  the  IIP  is  within  DCL  of  a  DISS  site  but  further  than  DCL  from  the  nearest  SSIES 
datum,  then  the  TEC  measurement  is  converted  into  an  equivalent  SSIES  density  measurement. 
First,  the  A <f0F2  and  A hmF2  corrections  determined  from  the  nearest  DISS  site  is  applied  to  the 
PRISM  profile  at  the  HP.  Then  the  topside  portion  of  the  profile  is  corrected  to  force  agreement 
with  the  TEC  measurement: 


(40) 


t(z)= 


Z, 

KFi+a{z-KFi\ 


z<hmF2 

z  >  KF2 


(41) 


where  a  is  a  scale  factor  to  be  determined  from  TEC. 
The  difference  in  TEC  is 


ATEC  =  f°°  rCT]dz-  f  n(*d)dz 

hmF,  O  JhmF,  O 


(42) 


or 


ATEC  =  — f  n(°ld)dz 
a  Kfi  ° 


(43) 


resulting  in 


a  =  ■ 


f  n^dz 

JhF,  O 


ATEC+f  n{°ld]dz 

Jh.F,  0 


(44) 
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If  we  identify  the  scale  factor  a  with  the  ratio  (z,  -hmF2)j  (zd  -hmF2),  then 


and 


=  KF2+{zd-KF. \)a 

(45a) 

N'0P  =”o'd)(Zl) 

(45b) 

II 

ft 

(45c) 

If  the  IIP  is  within  1000  km  of  both  a  DISS  station  and  a  DMSP  orbital  track,  the  TEC 
measurement  is  ignored.  In  a  future  version  of  PRISM,  in  which  the  plasmasphere  is  included, 
TEC  measurements  in  such  situations  will  be  used  to  constrain  the  H+  density. 


4.6  The  High  Latitude  Adjustment  Algorithm 

Due  to  the  complexity  of  the  high  latitude  ionosphere,  the  real  time  adjustment  algorithm 
differs  appreciably  from  the  low  and  midlatitude  algorithm.  Until  SSUSI  data  becomes 
available,  there  will  be  insufficient  data  to  adjust  the  parameterized  USU  model  in  the  way  that 
the  parameterized  low  and  midlatitude  models  can  be  adjusted.  Even  with  SSUSI  data,  it  is  not 
clear  that  the  midlatitude  algorithm  is  appropriate  for  the  high  latitude  regions. 

In  PRISM  1.2,  the  first  step  in  the  high  latitude  real  time  adjustment  process  is  the 
establishment  of  boundary  locations.  Three  boundaries  are  required:  (1)  the  equatorward  edge  of 
the  trough,  (2)  the  equatorward  edge  of  the  auroral  oval,  and  (3)  the  poleward  edge  of  the  auroral 
oval. 


The  equatorward  edge  of  the  trough  is  determined  from  SSIES  drift  meter  data  as  the 
point  where  the  measured  ion  drift  speed  departs  from  the  corotation  value.  The  trough  boundary 
as  a  function  of  magnetic  longitude  is  given  by  the  formula 


9,{<p)  =  O0(lt)  +  aexv 


M 


(46) 


where  cp  is  the  magnetic  local  time  (MLT,  hours),  and  It  is  a  "trough  index"  correlated  with  Kp. 
60{I, )  is  the  radius  of  the  trough  boundary  at  magnetic  local  midnight  and  is  given  by 


60(I,)  =  24.4°  +  2.12°/,  (47) 

The  second  term  represents  the  dayside  distortion  of  the  boundary,  which  would  otherwise  be  a 
circle  centered  on  the  magnetic  pole.  The  parameter  values  are 
a  =  -10.5° 
b  =  11.5  hr 
c  =  3.88  hr 
p  =  2.73 
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This  model  is  an  approximation  to  the  convection  boundaries  shown  in  Heppner  and  Maynard 
[1987]. 

The  trough  index,  /,,  is  determined  according  to  the  following  algorithm. 

1 .  If  there  is  no  ion  drift  data,  or  if  a  boundary  cannot  be  identified  in  the  data,  then  It~Kp. 

2.  If  a  single  boundary  crossing  is  identified  at  colatitude  6b  and  local  tim e<ph ,  then 


and 


90  =  6h-aexp 


j  &o  ~  24.4° 
'  2.12° 


(48) 

(49) 


3.  If  two  or  more  crossings  are  identified,  then  the  value  of  It  used  is  the  average  of  the 

values  determined  for  each  crossing. 

The  boundaries  of  the  auroral  oval  are  determined  from  electron  and  ion  precipitation 
data  from  the  SSJ/4  instrument.  Separate  boundaries  are  determined  for  electron  and  ion 
precipitation.  The  algorithm  for  determining  the  electron  and  ion  precipitation  boundaries  from 
the  SSJ/4  data  is  described  in  Appendix  D. 

The  next  step  in  the  high  latitude  adjustment  process  depends  on  the  amount  and  kind  of 
data  available.  The  decision  of  how  to  proceed  is  made  separately  for  the  E  layer  (NO+and  Oj) 
and  the  F  layer  (0+).  In  each  case,  two  choices  are  available: 

F  layer: 

1 .  Perform  a  simple  least  squares  adjustment  of  the  USU  0+  model. 

2.  Use  a  semi-empirical  foF2  model  (FMODEL)  to  adjust  the  USU  0+  profiles. 

E  layer: 

1.  Perform  a  simple  least  squares  adjustment  of  the  USU  NO+  and  Oj  models. 

2.  Use  a  fast,  first  principles,  E  layer  local  chemistry  model  (HLE). 


The  decision  matrix  used  in  PRISM  Table  5.  PRISM  High  Latitude  Decision  Matrix 


is  shown  in  Table  5. 

Data  Available? 

Model  Used  by  PRISM  | 

DISS 

SSIES 

SSJ/4 

F  layer 

E  layer 

In  experimenting  with  high 

Yes 

Yes 

Yes 

FMODEL 

HLE 

latitude  data,  we  found  that 

Yes 

Yes 

No 

FMODEL 

USU 

extrapolating  SSJ/4  data  taken  along 

Yes 

No 

Yes 

FMODEL 

HLE 

the  DMSP  orbital  track  to  points 

Yes 

No 

No 

USU 

USU 

well  away  from  the  orbital  track  was 

No 

Yes 

Yes 

USU 

USU 

very  risky.  We  found  no  suitable 

No 

Yes 

No 

USU 

USU 

model  of  the  instantaneous  auroral 

No 

No 

Yes 

USU 

USU 

precipitation  for  this  extrapolation. 

No 

No 

No 

USU 

USU 
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When  SSUSI  auroral  image  data  becomes  available,  this  limitation  will  be  removed  because 
much  less  extrapolation  will  be  required.  It  should  be  possible  to  use  HLE  whenever  timely 
SSUSI  images  are  available. 

FMODEL  is  a  semi-empirical  model  of  f0F2  based  on  a  combination  of  theory  and  data. 
It  is  divided  into  three  regions:  the  subauroral  trough,  the  auroral  oval,  and  the  polar  cap.  The 
f0F2  determined  by  least  squares  adjustment  of  the  model  parameters  is  used  to  scale  the  USU 

0+  profiles.  No  further  adjustment  of  the  profiles  is  performed. 

The  subauroral  trough  is  divided  into  two  local  time  regimes:  evening  (from  1200  to 
0000  MLT)  and  morning  (from  0000  to  1200  MLT).  Each  trough  (morning  and  evening)  has  a 
depth  parameter  that  specifies  the  difference  between  the  midlatitude  value  of  f0F2  at  the 
equatorward  edge  and  the  trough  minimum.  The  local  time  variation  of  the  trough  depth  is  fixed 
(not  part  of  the  least  squares  adjustment  process).  If  the  width  of  the  trough  is  less  than  3°,  the 
depth  is  reduced  in  proportion  to  the  width  so  that  when  the  width  vanishes  so  does  the  depth. 
The  thickness  of  the  (equatorward)  trough  wall  is  always  60%  of  the  total  width  of  the  trough. 
The  poleward  edge  of  the  trough  is  the  equatorward  edge  of  the  auroral  F-layer,  so  the  poleward 
"wall"  is  considered  to  be  part  of  the  auroral  region. 

At  fixed  magnetic  local  time,  the  auroral  F  layer  is  simply  a  cubic  polynomial  in 
magnetic  latitude: 


fM» = /«, + 4*-  a-)1  +*(*- <5°) 

where  fmax  =  f0F2(Xmax)  is  an  extremum,  and  A  and  B  are  chosen  so  that  f0F2  is  continuous 
across  the  boundaries  with  the  trough  and  polar  cap. 

The  background  polar  cap  f0F2  is  obtained  from  the  URSI  coefficients  using  an  internally 
derived  effective  sunspot  number,  the  value  of  which  is  determined  as  part  of  the  least  squares 
adjustment  process  based  on  ionosonde  data.  It  should  be  noted  that  this  part  of  the  model 
describes  only  the  background  polar  cap  ionosphere,  not  "polar  cap  patches"  or  "polar  cap  arcs." 


5.  VALIDATION 

AFRL  is  conducting  an  extensive  effort  to  validate  PRISM  using  some  of  the  many  new 
data  sources  that  have  become  available  in  recent  years.  The  role  of  CPI  in  this  effort  will  be 
documented  in  a  future  report. 
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6.  DISCUSSION 


Although  PRISM  represents  the  first  ionospheric  model  to  attempt  near  real  time 
assimilation  of  a  wide  variety  of  ionospheric  data,  it’s  development  began  a  decade  and  a  half 
ago.  Many  of  the  design  decisions  were  made  with  the  limitations  on  disk  storage,  computer 
memory,  and  CPU  processing  speeds  that  were  prevalent  at  the  time.  All  of  these  limitations 
have  been  lifted  and  those  design  decisions  should  be  reexamined.  For  example,  there  is  no 
longer  any  reason  not  to  increase  the  longitude  resolution  of  the  low  latitude  region  to  match  that 
of  the  midlatitude  region.  Nor  is  there  any  reason  not  to  replace  PRISM’s  seasonal  coefficient 
sets  with  monthly  coefficient  sets.  And  PRISM’s  operation  could  be  streamlined  by  using  a 
single  set  of  EOFs  for  all  latitudes,  seasons,  and  solar/geomagnetic  conditions.  This  would  make 
it  possible  to  interpolate  on  these  variables  in  coefficient  space  rather  than  configuration  space. 

Beyond  bringing  PRISM  up  to  date  with  modem  computer  capabilities,  there  are  a 
number  of  other  design  decisions  that  ought  to  be  revisited.  One  of  these  is  the  design  of  the 
Real  Time  Adjustment  (RTA)  algorithm.  While  it  has  the  virtue  of  reproducing  the  input  data 
(which  virtue  depends  on  the  integrity  and  accuracy  of  such  data),  it  is  not  capable  of  handling 
vertical  TEC  data  very  well,  and  slant  TEC  not  at  all.  The  use  of  a  Kalman  filter  or  Discrete 
Inverse  Theory  approach  ought  to  be  investigated  to  see  if  either  or  both  could  more  effectively 
facilitate  the  assimilation  of  slant  TEC  data. 

Another  design  feature  that  needs  to  be  reexamined  is  the  entire  High  Latitude  RTA. 
This  was  the  first  part  of  PRISM  to  be  developed  and  has  not  been  either  adequately  validated  or 
reexamined  since  that  time.  Kevin  Scro  of  AFSMC,  who  is  responsible  for  transitioning 
software  to  AFWA,  has  expressed  legitimate  concerns  over  the  fact  that  the  algorithm  changes 
discontinuously  depending  on  the  kind  of  data  available  to  PRISM.  The  entire  high  latitude 
algorithm  needs  to  be  reexamined  and  reformulated  to  take  into  account  the  capability  of 
mapping  the  auroral  oval  boundaries  over  much  of  the  high  latitude  region  from  SSUSI  images 
that  will  soon  be  forthcoming  from  DMSP.  Other  enhancements  that  should  be  considered  are 
(a)  the  addition  of  a  good  plasmasphere  model,  whether  empirical  or  theoretical  and  (b)  the 
inclusion  of  low  latitude  vertical  drift  (especially  the  pre-reversal  enhancement)  among  the 
geophysical  parameters. 
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Appendix  A.  Empirical  Orthonormal  Functions 

This  treatment  of  empirical  orthogonal  functions  (EOF's)  is  based  on  the  Appendix  of 
Secan  and  Tascione  [1984],  which  was  based  on  Lorenz  [1956],  Kutzbach  [1967],  and  Davis 
[1976].  See  also  Peixota  and  Oort  [1991].  The  reader  is  referred  to  these  references  for 
mathematical  proofs  of  the  assertions  made  below.  In  the  following  discussion,  we  use  the 
notation  given  in  Table  3  of  the  main  text. 

A  database  consists  of  altitude  profiles  at  certain  longitudes,  certain  latitudes,  and  certain 
Universal  Times.  (See  Tables  1-4  of  the  main  text.)  Let  S  be  the  number  of  altitude  profiles  in  a 
database,  and  let  /  be  the  number  of  points  in  each  altitude  profile.  We  would  like  to  represent 
each  altitude  profile  of  the  quantity  ¥  (e.g.,  0+ concentration)  as  an  expansion  in  orthogonal 
functions,  gm(z,): 


^M  =  Ea-^.(zi)  +  r.(zi)»  s  =  l..S,  /  =  !.../ 


m=\ 


(Al) 


where  rs{zt)  is  the  residual,  and  the  coefficients  asm  are  calculated  from 

=  !*,(*,)&(*,)  (A2) 

1-1 

In  principle,  any  orthogonal  set  of  functions  may  be  used.  However,  the  references  cited  above 
provide  an  algorithm  for  finding  the  set  which  minimizes  the  RMS  error  for  a  given  number  of 
terms,  M<  /.  We  summarize  the  algorithm  here. 

First  define  the  /  x  /  covariance  matrix  C  with  elements 


c6=jZ'!'.M'p.(4 

O  j=l 


i,j  =  1,2,.../ 


(A3) 


Now  consider  the  eigenvalue/eigenvector  problem  Cep  =  cpL  or 

Z  c„<pA  =  ■  VA  <A4> 

J= l  M 

where  cp  =  {#>,., }  is  the  matrix  of  eigenvectors  of  C  =  {dZ},  and  L  =  is  a  diagonal  matrix 

whose  elements  are  the  corresponding  eigenvalues.  (The  A*  column  of  cp  is  the  eigenvector 
corresponding  to  the  eigenvalue,  Xk.)  By  convention,  the  eigenvectors  and  eigenvalues  are 
ordered  so  that  A,  >  X2  >...>  Xr  Because  C  is  a  real  symmetric  matrix,  eigenvectors 
corresponding  to  unique  eigenvalues  are  guaranteed  to  be  orthogonal  [See,  e.g.,  Hildebrand, 
1965].  Because  of  the  origin  of  the  matrix  C,  it  is  unlikely  that  any  of  its  eigenvalues  will  be 
degenerate,  so  we  may  assume  that  cp  is  an  orthogonal  set.  According  to  Secan  and  Tascione 
[1984]  and  references  therein,  the  set  of  orthogonal  functions  that  minimizes  the  RMS  error  for 
M  terms  is  just  the  first  M  eigenvectors: 

gm(z,.)  =  ^„  i  =  l,2,.../;  m  =  1,2,... M  (A5) 


These  are  the  Empirical  Orthonormal  Functions  (EOFs). 

As  a  practical  matter,  we  have  found  that  the  number  of  EOFs  needed  to  provide  a 
reasonably  good  representation  for  all  the  profiles  is  about  // 6,  as  illustrated  in  Table  4  in  the 
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main  text.  The  only  exception  is  the  low  and  midlatitude  E  layer  (NO+  and  Oj),  probably 
because  these  databases  covered  both  hemispheres  simultaneously.  We  have  also  found  that 
substantial  improvement  in  representation  does  not  occur  until  the  number  of  EOF's  is  about  III. 
Furthermore,  the  EOF's  derived  for  one  database  were  inadequate  for  any  other  database,  and  the 
EOF's  simultaneously  derived  from  several  databases  produce  noticeably  poorer  representations 
than  those  derived  for  each  database  individually.  Consequently,  we  have  derived  separate  EOF 
sets  for  each  database. 

The  first  nine  EOF's  derived  from  the  low  latitude  F  region  (0+)  database  for  the  US 
longitude  sector,  the  December  solstice,  and  moderate  solar  activity,  are  shown  in  Figure  Al. 
The  first  EOF  always  has  the  least  structure,  and  successive  EOF's  become  progressively  more 
structured.  Although  differing  in  detail,  the  EOF’s  for  the  other  databases  are  qualitatively 
similar. 


Low  Latitude  0* 


Figure  Al.  The  Empirical  Orthonormal  Functions  (EOF's)  for  low  latitude  0+  derived  from  the  LOWLAT 
output  databases  for  the  USA  longitude  sector,  December  solstice,  moderate  magnetic  activity, 
and  moderate  solar  activity.  Only  the  first  nine  EOF’s  are  plotted  because  these  are  the  ones  used 
in  PRISM. 
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Appendix  B.  Orthogonal  Polynomials  of  Discrete  Variables 

Because  the  databases  to  which  we  desire  analytic  approximations  have  discrete  latitude 
grids,  we  preferred  to  use  polynomials  whose  orthogonality  is  defined  in  terms  of  that  grid,  rather 
than  in  terms  of  integrals  over  the  interval.  The  algorithm  for  generating  orthogonal  polynomials 
on  a  specified  grid  is  given  by  Beckmann  [1973].  Let  us  denote  the  desired  polynomials  by  un  (X) 
and  define  m.,(X)  =  0  and  u0(K)  =  1.  Note  that  the  polynomials  are  continuous  functions  of  the 
continuous  variable  A  even  though  their  orthogonality  is  defined  in  terms  of  the  discrete  grid 
| Aj,j  =  1,2,. ..jj.  The  recursion  relation  for  the  polynomials  is 

(Bi) 

where  the  norms  hn  are  given  by 

K l02) 

7=1 

and  the  recursion  constants  Bn  are  given  by 

<B3> 

nn  7=1 

The  polynomials  generated  by  this  algorithm  may  be  used  to  represent  the  latitude 


variations  of  the  Fourier  coefficients  amp  and  bmp  (see  main  text): 

N 

ar,v{xJ’T)  =  Za*v(r/K(A;) 

n=\ 

(B4) 

Kp^jSl )  =  S  Pmnp iTl  )Un  (^J ) 
n- 1 

(B5) 

where 

^(r/)  =  p-Za^(A  y»r/K(A/) 

nn  j-l 

(B6) 

Pmnp^i)  =  ^7lLK{XPZ^n{^j) 

nn  M 

(B7) 
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Appendix  C.  Auroral  boundary  determination 


There  are  four  auroral  precipitation  boundaries  in  each  hemisphere:  The  equatorward  and 
poleward  electron  boundaries  and  the  equatorward  and  poleward  ion  boundaries.  Since  there  are 
always  two  DMSP  satellites  in  orbit,  there  are  potentially  four  crossings  of  each  boundary.  The 
specification  of  a  given  boundary  depends  on  the  number  of  crossings  detected  in  the  data.  The 
boundary  detection  is  performed  as  follows. 

1.  If  SSJ/4  data  are  available,  then  divide  each  orbit's  worth  of  data  into  four  segments,  each 
segment  extending  from  the  most  equatorward  point  to  the  most  poleward  point. 

2.  Each  segment  is  searched  from  equator  toward  the  pole  until  the  electron  energy  flux  exceeds 
0.25  erg  cm'2  s'1  (for  the  equatorward  electron  boundary)  and  the  ion  energy  flux  exceeds 
0.1  erg  cm'2  s'1  (for  the  equatorward  ion  boundary). 

3.  Search  each  segment  from  pole  toward  the  equator  until  the  same  thresholds  are  exceeded  to 
establish  the  poleward  boundaries. 

The  boundaries  are  always  assumed  to  be  circular.  The  specification  of  each  boundary  circle  (in 
terms  of  its  center  and  radius)  depends  on  die  number  of  crossings  detected. 

1.  If  four  crossings  are  found,  there  are  four  possible  combinations  of  three  points.  Each  of  the 
four  possible  combinations  defines  a  circle,  which  may  be  described  by  its  center  and  radius. 
The  center  coordinates  and  radii  of  the  four  circles  are  averaged  to  determine  the  boundary 
used  in  PRISM. 

5.  If  only  three  crossings  are  found,  the  boundary  is  uniquely  defined  by  the  circle  passing 
through  the  three  points. 

6.  If  only  two  crossings  are  found,  the  boundary  is  taken  to  be  the  circle  that  passes  through 
both  points  and  has  a  radius  equal  to  the  mean  of  the  colatitudes  of  the  two  points. 

If  only  one  crossing  of  a  boundary  is  found,  or  if  no  crossings  are  found  (or  SSJ/4  data  is 
missing),  the  boundaries  are  determined  as  follows.  Equatorward  boundaries  are  based  on  an 
analytic  representation  of  the  Gussenhoven  et  al.  [1983]  boundaiy.  The  boundary  itself  is  a  circle 
whose  center  is  displaced  from  the  magnetic  pole.  The  radius  (in  degrees)  of  the  circle  is 


0,  (/>)  =  20.9°  + 1.7° />  (Cl) 

where  Px  is  a  "precipitation  index".  The  center  of  the  circle  is  located  at 

JLct  (/> )  =  87.3°-0.267°  P,  (C2) 

tpcl  (/> )  =  39.5°-125°  />  +  0.076°  />2  (C3) 

where  A,cI  is  magnetic  latitude  and  (pc]  is  magnetic  local  time  (in  degrees). 
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Poleward  boundaries  have  almost  the  same  form,  except  that  they  are  parameterized  in 
terms  of  a  separate  "precipitation  index"  P2. 


02{P2)  =  13.4°  + 1.7°  P2 

(C4) 

J89.2°+0.267°P2,  P2<  3 

Xcl  > =  |90.8°-0.267°  P2,  P2>  3 

(C5) 

f^  +  1800,  P2<  3 

(C6) 

If  there  is  no  data,  or  if  no  boundaries  are  detected  in  the  data,  then/^  =  P2=Kp.  However,  if  a 
single  crossing  is  detected,  then  the  precipitation  index  is  chosen  so  that  the  above  boundary 
matches  the  crossing. 
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Appendix  D.  PRISM  Input/Output  File  Specifications 

PRISM  1 .7e  External  Input/Output  File  Specifications 

(for  files  created  by  or  under  the  control  of  the  User) 
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This  record  is  present  only  if  FNORM  (record  7)  is  0. 

This  record  is  present  only  if  STjOUT  (record  5)  is  0, 1,  or  2. 
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The  Data  Element  Names  do  not  necessarily  correspond  to  the  FORTRAN  variable  names  in  the  source  code.  This  is  because  the  FORTRAN 
names  are  not  always  descriptive  of  the  actual  contents  of  the  data  element  (for  historical  reasons). 

Because  PRISM  uses  FORTRAN  list  directed  input  for  all  numeric  input,  the  numeric  format  is  flexible.  The  formats  indicated  in  the  table  are 
suggested  formats  only.  Plus  (+)  signs  are  always  optional.  The  variables  may  be  separated  by  (1)  one  or  more  spaces,  (2)  a  comma  (,),  or  (3)  a 
comma  and  one  or  more  spaces. 
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All  records  have  the  same  format.  PRISM  recognizes  the  end-of-file  (eof)  mark  as  an  “end-of-data”  indicator. 
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Changed  maximum  allowed  number  of  accepted  SSIES  IN  SITU  PLASMA  records  in  Note  6;  added  explanation 
for  that  number  to  Note  6; 

Added  explanation  of  maximum  allowed  number  of  accepted  SSJ/4  records  to  Note  10) 
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Appendix  F.  Changes  Memos  for  Versions  1.6  through  1.7e 


The  following  memoranda  provide  detailed  descriptions  of  the  changes  made  to  each  version  of 
PRISM  to  produce  the  next  version. 
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MEMORANDUM 


DATE:  12-August-1995 

TO:  RobDaniell 

FROM:  Lincoln  Brown 

RE:  Changes  to  PRISM  1.5  for  PRISM  1.5a 


The  changes  to  PRISM  1.5  for  PRISM  1.5a  focus  on  a  change  to  the  Station  Output  File  Specification  requested  by 

Hughes  STX,  the  correction  of  reported  problems  regarding  ingestion  and  use  of  DMSP  and  TISS  data,  and 

miscellaneous  bug  fixes.  The  changes  are  summarized  as  follows: 

I.  The  third  header  line  at  the  beginning  of  the  station  output  file  has  been  removed  and  a  header  line  added 
above  each  critical  parameters  data  record  for  station  output  type  IOUTS=0.  This  was  requested  by  Hughes  STX 
and  requires  a  change  to  the  Station  Output  File  Specification. 

n.  The  ingestion  of  DMSP  data  has  been  corrected  and  strengthened  based  on  problems  reported  by  Hughes 
STX.  PRISM  can  now  handle  up  to  24  sets  of  each  kind  of  DMSP  data  in  a  DMSP  data  file,  and  it  can  handle  up 
to  8  DMSP  data  files. 

III.  The  ingestion  of  TEC  data  has  been  corrected  based  on  problems  reported  by  George  Bom’s  group. 

IV.  Under  certain  conditions,  several  uninitialized  variables  in  common  block  INDIRECT  were  referenced. 
They  are  now  initialized  in  BLOCK  DATA  ITR4. 

V.  The  internal  option  of  producing  0+  density  instead  of  electron  density  on  output  has  been  removed  because 
it  was  not  correctly  implemented. 

VI.  Several  subroutines  that  are  no  longer  used  have  been  removed. 

The  table  below  describes  the  changes  that  I  made  to  PRISM  1.5  to  produce  PRISM  1.5a. 
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Module 

Description  of  Changes  (Continued) 

OUTPUT.FOR 

Subroutine  WR_ST_DATA 

Removed  the  third  header  line  at  the  beginning  of  the  station  output  file  for  station  output  type 
IOUTS=0. 

Added  a  header  line  above  the  critical  parameters  data  record  line  in  the  station  output  file  for  station 
output  type  IOUTS=0. 

PRECIP.INC 

n/a 

Added  PARAMETER  MAXORB. 

Changed  the  last  dimension  of  common  block  PRECP1  arrays  LAT,  MLT,  and  ERG  from  2  to 
MAXORB  to  support  more  DMSP  SSJ/4  data. 

Changed  the  last  dimension  of  common  block  PRECP1  arrays  MLATR,  MLTR,  and  VALTR  from  2 
to  MAXORB  to  support  more  DMSP  Ion  Drift  data. 

PRISM.  FOR 

Block  Data  ITR4 

Added  initialization  of  common  block  INDIRECT  variables  HDFE,  HDHE,  HDF2,  and  HDHF2. 

Block  Data  PREC 

Complete  rewrite  for  more  flexible  initialization  of  arrays  in  common  block  PRECP1 . 

ES13S323S!3HHIHHI^IH 

Updated  the  version  number  and  version  date. 

j  Ends  | 

MEM  OR- AN*  D-  EM 


DATE: 

29-September-l  995 

TO: 

Rob  Daniell 

FROM: 

Lincoln  Brown 

RE: 

Changes  to  PRISM  1.5a  for  PRISM  1.6 

The  changes  to  PRISM  1 .5a  for  PRISM  1 .6  focus  on  the  resolution  of  problem  report  PRF-PRISM 1 9  from  Bob 
Prochaska  at  Hughes  STX  and  final  enhancements  for  the  last  planned  delivery  of  PRISM  1 .  The  changes  are 
summarized  as  follows: 

I.  Problem  report  PRF-PRISM1 9  has  been  resolved.  PRISM  uses  the  plasma  temperature  to  calculate  a 
topside  scale  height  for  the  midlatitude  electron  density,  where  the  plasma  temperature  is  the  sum  of  the  electron 
and  ion  temperatures  from  a  DMSP  SSIES  In  Situ  Plasma  data  record.  The  topside  scale  height  varies  linearly 
with  the  plasma  temperature.  Missing  electron  and  ion  temperature  data  in  the  DMSP  SSIES  In  Situ  Plasma  data 
record  is  flagged  by  zero  values.  Previously,  however,  if  one  of  the  temperatures  was  zero  but  the  other 
temperature  was  nonzero,  PRISM  continued  to  calculate  and  use  a  plasma  temperature,  resulting  in  a  topside 
scale  height  roughly  half  as  large  as  it  should  have  been.  The  erroneously  small  scale  height  drove  the  topside 
midlatitude  electron  density  to  fall  off  rapidly,  resulting  in  very  a  small  electron  density  at  the  top  of  the  altitude 
grid.  Now,  PRISM  ignores  the  temperature  data  in  a  DMSP  SSEES  In  Situ  Plasma  data  record  and  does  not 
calculate  a  topside  scale  height  for  die  data  record  if  either  the  electron  or  ion  temperature  is  zero. 

II.  PRISM  now  checks  Ionosonde  (DISS)  data  for  nonphysical  values.  The  following  tests  have  been  added  to 
the  Ionosonde  data  ingestion: 

A.  The  observed  critical  frequency  of  the  F2  layer  (foF2)  in  the  Ionosonde  data  record  must  be  in  the  range 
0  <foF2  <  28.4  MHz.  The  range  maximum  corresponds  to  a  peak  F2  layer  density  (nJS)  of  1 07  cm'3  using 
the  approximation  KmF2(cm'3)=l  ,24x  104-/oF2(MHz)2.  If  the/oF2  parameter  in  the  Ionosonde  data  record  is 
outside  this  range,  then  the  foF2  parameter  is  ignored. 

B.  The  observed  height  of  the  F2  layer  (h„F2)  in  the  Ionosonde  data  record  must  be  in  the  range 

200  <  hJF2  <  1000  km.  If  the  A„F2  parameter  in  the  Ionosonde  data  record  is  outside  this  range,  then  the 
hmF2  parameter  is  ignored. 

C.  The  observed  critical  frequency  of  the  E  layer  (fJT)  in  the  Ionosonde  data  record  must  be  in  the  range 

0  <foE  <  28.4  MHz.  The  range  maximum  corresponds  to  a  peak  E  layer  density  («,„£)  of  107  cm'3  using  the 
approximation  n„E( cm"3)=l  .24x1 04  /^(MHz)2.  If  the fE  parameter  in  the  Ionosonde  data  record  is  outside 
this  range,  then  the  /^/s  parameter  is  ignored. 

D.  The  observed  height  of  the  E  layer  (hJEJ)  in  the  Ionosonde  data  record  must  be  in  the  range  90  <  hJL  < 

1 50  km.  If  the  hmE  parameter  in  the  Ionosonde  data  record  is  outside  this  range,  then  the  hmE  parameter  is 
ignored. 

TTT  PRISM  now  checks  TEC  (IMS)  data  for  nonphysical  values.  The  vertical  equivalent  TEC  in  the  TEC  data 
record  must  be  in  the  range  0  <  TEC<  250  TEC  Units  (1  TEC  Unit=1012  cm‘2=  10 16  m'2)  .  If  the  vertical 
equivalent  TEC  in  the  TEC  data  record  is  outside  this  range,  then  the  TEC  data  record  is  ignored. 
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IV.  PRISM  now  checks  DMSP  SSIES  In  Situ  Plasma  data  for  nonphysical  values.  The  following  tests  have 
been  added  to  the  DMSP  SSIES  In  Situ  Plasma  data  ingestion: 

A.  The  in  situ  electron  density  ( ne )  in  the  DMSP  SSIES  In  Situ  Plasma  data  record  must  be  in  the  range 
103  <ne<  107  cm'3.  If  the  ne  parameter  in  the  DMSP  SSIES  In  Situ  Plasma  data  record  is  outside  this  range, 
then  the  entire  data  record  is  ignored. 

B.  The  in  situ  electron  temperature  (Te)  in  the  DMSP  SSIES  In  Situ  Plasma  data  record  must  be  in  the 
range  0  <Te<  6000  °K.  The  in  situ  ion  temperature  (7})  in  the  DMSP  SSIES  In  Situ  Plasma  data  record 
must  be  in  the  range  0  <Tt<  6000  °K.  If  either  the  electron  or  ion  temperature  in  the  DMSP  SSIES  In  Situ 
Plasma  data  record  is  outside  its  range,  then  both  temperatures  are  ignored. 

V.  Phantom  Ionosonde  data  has  been  added  to  PRISM.  In  the  absence  of  any  valid  real-time  Ionosonde  and 
TEC  data  in  30°x30°  corrected  geomagnetic  latitude/longitude  zones  between  -60°  and  +60°  corrected 
geomagnetic  latitude,  PRISM  now  adds  phantom  Ionosonde  data  centered  in  the  zones  to  reduce  the  global  effect 
of  sparse  real-time  Ionosonde  and  TEC  data.  The  phantom  Ionosonde  data  is  generated  using  the  parameterized 
models  in  PRISM  for  the  UT  of  the  run,  and  is  included  in  the  set  of  real-time  data  that  PRISM  uses  for  its  real¬ 
time  adjustment  of  the  parameterized  models. 

VI.  Phantom  In  Situ  Plasma  data  has  been  added  to  PRISM.  In  the  absence  of  any  real-time  SSIES  In  Situ 
Plasma  data  in  30ox30°  corrected  geomagnetic  latitude/longitude  zones  between  -60°  and  +60°  corrected 
geomagnetic  latitude,  PRISM  now  adds  phantom  In  Situ  Plasma  data  (density  only  since  PRISM  does  not  contain 
electron  and  ion  temperature  models)  centered  in  the  zones  to  reduce  the  global  effect  of  sparse  real-time  In  Situ 
Plasma  data.  The  phantom  In  Situ  Plasma  density  data  is  generated  using  the  parameterized  models  in  PRISM 
for  the  UT  of  the  run  and  DMSP  altitude  (840  km),  and  is  included  in  the  set  of  real-time  data  that  PRISM  uses 
for  its  real-time  adjustment  of  the  parameterized  models. 

VII.  Commented-out  coding  has  been  removed. 

VIII.  Unused  routines  have  been  removed. 

IX.  Unused  variables  have  been  removed. 

The  table  below  describes  the  changes  that  I  made  to  PRISM  1 .5a  to  produce  PRISM  1 .6. 


Module 

Program  Unit 

Description  of  Changes  (Begins) 

FMODEL.FOR 

Subroutine  FMODEL 

Removed  local  variables  TFOF2,  THMF2,  CFOF2,  POFOF2,  CFOE,  and  CHME  since  they  are  not 
used. 

Subroutine  TRO_DEP 

Removed  local  variables  FOF2,  HMF2,  FOF2ML,  and  HMF2ML  since  they  are  not  used. 

GETDAT.FOR 

Subroutine  DOJDIR 

Removed  argument  NDUMDAT(I,2)  from  the  first  call  to  routine  DO_JONO  since  it  is  not  used  by 
that  routine. 

Removed  argument  NDUM1  from  the  second  call  to  routine  DO _JONO  since  it  is  not  used  by  that 
routine. 

Removed  local  variable  NDUM1  since  it  is  no  longer  used. 

Subroutine  DO.JONO 

Expanded  the  range  validation  of  ingested  E  and  Fi  layer  Ionosonde  data. 

Removed  commented-out  code  for  defunct  BOTTOMSIDE  data  type. 

Removed  PARAMETER  AFPE  since  it  is  no  longer  used. 

Removed  output  parameter  NBREC  since  it  is  no  longer  used. 

Subroutine  DO_TEC 

Expanded  the  range  validation  of  ingested  TEC  data. 

Subroutine  DCLDMSP 

Expanded  the  range  validation  of  ingested  DMSP  SSIES  In  Situ  Plasma  data.  This  resolves  problem 
report  PRF-PRISM19  from  Bob  Prochaska  at  Hughes  STX. 

Removed  PARAMETER  TDIFMAX  since  it  is  no  longer  used. 

Subroutine  IES_DATA 

Removed  since  it  is  not  used. 

ICLUTIL.FOR 

Subroutine  CLOSFL 

Removed  since  it  is  not  used. 

|  Continues  | 
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Module 


MATHJJTI.FOR 


Program  Unit 


Subroutine  SORT7 


Subroutine  SWPARR 


Subroutine  INDEXX 


Subroutine  GETGAM 


Subroutine  MIDLAT 


Subroutine  PHIONO 


Subroutine  DST1 


Subroutine  PHIES 


Subroutine  INS  TRACK 


n/a 


Program  PRISM 


Description  of  Changes  ( Continue 


MIDLAT.FOR 


PHANTOM.  FOR 


PRISM. FOR 


READ  DBA.FOR 


RTA.FOR 


TIMELIB.FOR 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


READAWS 


RTA 


BOTTOM 


STRDEL 


STRFLL 


STRINS 


STRNPO 


STRRPL 


STRSHC 


STRSHI 


STRTRM 


STRUCA 


CONCAT 


DELETE 


FILL 


INSERT 


LOCASE 


LENGTH 


NTHPOS 


REPLAC 


SHIFTC 


SHIFT 


TRIM 


UPCASE 


TIMMDT 


DELTAT 


GETDOY 


LEAPYR 


GETMDM 


MODATE 


Removed  since  it  is  not  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used. 


Removed  commented-out  PRINT  statements. 

Removed  FORMAT  statement  91000  since  it  is  no  longer  used. 

Removed  local  variables  AF  and  FIRST  since  they  are  no  longer  used.  _ _ _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  not  used. 


Removed  since  it  is  no  longer  used. 


New  module  containing  routine  PHANTM. 


Added  calls  to  routines  INIT  and  PHANTM  after  the  conversion  of  the  nominal  UT  from  hours  to 
seconds. 

Undated  the  version  number  and  version  date. 


Commented  out  references  to  local  variable  KF000  since  it  is  not  needed.  _ 


Removed  commented-out  code  for  bottomside  fittm 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  no  longer  used.  _ _ 


Removed  since  it  is  no  longer  used.  _ _ 


Removed  since  it  is  no  longer  used.  _ _ 


Removed  since  it  is  no  longer  used.  _ _ 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used.  _ _ _ 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  not  used. 


Removed  since  it  is  not  used. 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used. 


Removed  smce  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used. _ 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ 


Removed  since  it  is  not  used.  _ _ 


Removed  since  it  is  not  used.  _ 


Ends  _ _ 


MEMORANDUM 


DATE: 

12-February-1996 

TO: 

Rob  Daniell 

FROM: 

Lincoln  Brown 

RE: 

Changes  to  PRISM  1. 6  for  PRISM  1. 6a 

The  changes  to  PRISM  1 .6  for  PRISM  1 .6a  focus  on  the  replacement  of  the  LLF  parameterized  model  coefficients, 
the  resolution  of  PRF  PRISM-24,  and  a  change  to  the  Gridded  Output  Specification.  The  changes  are  summarized  as 
follows: 

I.  The  LLF  parameterized  model  coefficients  have  been  replaced,  resulting  in  several  improvements  in  the 
low-latitude  0+  density  representation  in  PRISM: 

A.  The  solar  activity  dependence  is  now  correctly  represented.  Previously,  due  to  bugs  in  the  LOWLAT 
theoretical  model  used  to  generate  the  coefficients,  and  misuse  of  one  of  the  LOWLAT  input  files,  no 
variation  in  the  0+  density  due  to  solar  activity  level  was  present. 

B.  The  peak  of  the  dayside  0+  density  is  more  in  line  with  expected  dayside  densities.  It  is  believed  that 
the  dayside  electron  temperature  model  used  in  LOWLAT  resulted  in  dayside  0+  densities  that  were  too 
large.  The  upper  limit  of  the  dayside  electron  temperature  has  been  increased  from  2500°K  to  3000°K  to 
maintain  a  smaller  dayside  0+  density. 

C.  The  bottomside  0+  density  has  been  improved.  Previously,  spline  interpolation  of  rapidly  decreasing  0+ 
density  on  a  sparse  bottomside  altitude  grid  could  result  in  a  false  and  substantial  peak  near  the  bottom  of  the 
LLF  altitude  grid  (160  km).  The  interpolation  algorithm  has  been  modified  to  eliminate  this  problem. 

II  PRF  prisM-24  has  been  resolved.  DMSP  SSIES  Ion  Drift  and  SSJ/4  data  sets  are  now  no  longer  allowed 
to  overrun  allocated  space  in  PRISM.  Previously,  if  the  maximum  allowed  number  of  data  points  was  accepted 
by  PRISM  before  the  end  of  the  dataset  was  reached,  an  array-out-of-bounds  error  occurred.  A  slight  change  in 
logic  has  resolved  this  problem.  This  problem  was  reported  by  Bob  Prochaska  at  Hughes  STX. 

III.  The  coordinate  system  flag  IGM  in  the  Gridded  Output  Specification  is  now  consistent  with  the  PRISM 
source  code.  Previously,  the  convention  given  in  the  Gridded  Output  Specification  was  opposite  that  used  in 
PRISM.  This  problem  was  reported  by  Vince  Eccles  at  Space  Environment  Corp.,  and  its  resolution  was 
approved  by  Bob  Prochaska  at  Hughes  STX. 

The  table  below  describes  the  changes  that  I  made  to  PRISM  1 .6  to  produce  PRISM  1 .6a. 


Module 

Description  of  Changes 

GETDAT.FOR 

Subroutine  DOJDMSP 

Changed  the  line  "IF(NSSIES  .LE.  MAXIES)  GOTO  150"  to  "IF(NSSIES  .LT.  MAXIES)  GOTO 

1 50"  to  avoid  an  array-out-of-bounds  error  when  MAXIES  data  points  in  a  SSIES  Ion  Drift  data  set 
have  been  accepted  before  the  end  of  the  data  set  has  been  reached. 

Changed  the  line  "IF(NSSJ4  .LE.  MAXJ4)  GOTO  250"  to  "IF(NSSJ4  .LT.  MAXJ4)  GOTO  250"  to 
avoid  an  array-out-of-bounds  error  when  MAXJ4  data  points  in  a  SSJ/4  data  set  have  been  accepted 
before  the  end  of  the  data  set  has  been  reached. 

PRISM.  FOR 

Program  PRISM 

Updated  the  version  number  and  version  date. 
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M-EM-ORANDU-M 


DATE: 

19-August-1996 

TO: 

Rob  Daniell 

FROM: 

Lincoln  Brown 

RE: 

Changes  to  PRISM  1.6a  for  PRISM  1.6b 

The  changes  to  PRISM  1.6a  for  PRISM  1.6b  focus  on  improving  the  midlatitude  real-time  adjustment.  The  changes 

are  summarized  as  follows: 

I.  Several  improvements  have  been  made  to  the  midlatitude  real-time  adjustment: 

A.  Real-time  and  phantom  data  that  agree  exactly  with  the  parameterized  models  are  no  longer  prevented 
from  influencing  the  midlatitude  real-time  adjustment.  Previously,  data  agreeing  exactly  with  the 
parameterized  models  was  ignored  by  the  midlatitude  real-time  adjustment.  This  problem  was  reported  by 
David  Coxwell  at  Air  Force  Institute  of  Technology. 

B.  The  actual  magnetic  latitude  of  real-time  and  phantom  data  is  now  used  in  the  distance-weighting  for  the 
midlatitude  real-time  adjustment.  Previously,  the  magnetic  latitude  of  data  used  in  the  distance-weighting  was 
normalized  to  the  range  [-90,90]  degrees  based  on  the  magnetic  latitude  of  the  equatorward  trough  boundary. 
The  normalization  introduced  error  into  the  midlatitude  correction  field  by  artificially  extending  the 
latitudinal  distance  between  the  data  and  the  location  of  interest,  resulting  in  an  elliptical  rather  than  a  circular 
weighting  function. 

C.  The  actual  magnetic  longitude  of  real-time  and  phantom  data  is  now  used  in  the  distance-weighting  for 
the  midlatitude  real-time  adjustment.  Previously,  the  magnetic  local  time  of  the  data  converted  to  hour  angle 
was  used  in  the  distance-weighting.  The  use  of  the  magnetic  local  time  instead  of  magnetic  longitude 
introduced  as  much  as  a  30  degree  (2  hour)  error  in  magnetic  longitude  in  the  midlatitude  correction  field, 
depending  on  the  time  difference  between  the  data  and  the  nominal  UT. 

D.  The  peak  density  corrections  of  the  midlatitude  real-time  adjustment  have  been  improved  by  using  total 
ion  densities  at  hnE  and  hmF2  instead  of  the  molecular  ion  density  at  l^E  and  the  0+  density  at  hJF^  This 
removes  the  assumptions  that  0+  is  negligible  at  hmE  compared  to  the  molecular  ions  and  that  molecular  ions 
are  negligible  at  hniF2  compared  to  0+. 

E.  The  density  scaling  correction  of  the  midlatitude  real-time  adjustment  has  been  improved  by  applying 
the  E-layer  and  F2-layer  density  scalings  to  both  the  molecular  ion  density  profile  and  the  0+  density  profile, 
resulting  in  a  self-consistent  correction.  This  removes  the  assumptions  that  0+  is  negligible  at  l^E  compared 
to  the  molecular  ions  and  that  molecular  ions  are  negligible  at  hmF2  compared  to  0+. 

II.  Several  minor  changes  have  been  made  for  compatibility  with  Microsoft  FORTRAN.  They  do  not  impact 

the  results. 

m  Commented-out  code  has  been  removed  from  several  routines. 

IV.  Unused  FORMAT  statements  have  been  removed  from  several  routines. 

V.  Several  typographical  errors  in  comments  have  been  corrected. 


The  table  below  describes  the  changes  that  I  made  to  PRISM  1.6a  to  produce  PRISM  1.6b. 


Module 

Program  Unit 

Description  of  Changes 

- 

MIDLAT.FOR 

Subroutine  MIDLAT 

Removed  conditionals  so  that  data  agreeing  exactly  with  the  model  is  not  ignored  in  the  midlatitude 
real-time  adjustment. 

Magnetic  local  time  stored  in  local  arrays  IMLT  and  JMLT  is  no  longer  converted  to  hour  angle. 

Added  local  arrays  IMLON  and  JMLON. 

Magnetic  longitude  is  now  stored  in  local  arrays  IMLON  and  JMLON. 

Local  array  JMLON  is  now  passed  to  routine  PER_ARR. 

Corrected  typo  in  comment  by  changing  U,HME  =  8”  to  “,HME  =  6”. 

Subroutine  PER_,ARR 

Added  array  JMLON  to  argument  list. 

Magnetic  local  time  stored  in  argument  array  JMLT  is  no  longer  expressed  as  hour  angle. 

Magnetic  latitude  stored  in  argument  array  JMLAT  is  no  longer  normalized  to  the  range  [-90,90] 
degrees  by  the  magnetic  latitude  of  the  equatorward  trough  boundary. 

Magnetic  longitude  is  now  stored  in  argument  array  LON  instead  of  magnetic  local  time. 

OUTPUT.FOR 

Subroutine  WRUE_DATA 

Changed  “1  PEI  1.4”  to  “,1  PEI  1.4”  in  FORMAT  statement  8000  to  avoid  compiler  error  under 
Microsoft  Fortran  PowerStation. 

Removed  commented-out  code. 

Removed  unused  FORMAT  statements. 

Subroutine  W_ST_DATA 

Removed  commented-out  code. 

Removed  unused  FORMAT  statements,  including  FORMAT  statement  8000,  which  caused  a 
compiler  error  under  Microsoft  Fortran  PowerStation. 

PARAM.FOR 

Subroutine  PARAM 

Magnetic  longitude  is  now  passed  to  routine  RTA  instead  of  magnetic  local  time  by  changing 
argument  MLT  to  MLON  in  the  call  to  routine  RTA. 

Removed  commented-out  code. 

PHANTOM.  FOR 

Subroutine  PHANTM 

Corrected  typos  in  comments  by  changing  “to  added”  to  “to  be  added”. 

PRISM.  FOR 

Block  Data  INTRAT1 

Character  strings  in  array  LABEL  longer  than  the  allocated  length  have  been  truncated. 

The  DATA  statement  defining  the  LABEL  array  has  been  reformatted  to  eliminate  line  continuations 
in  the  middle  of  character  strings. 

Program  PRISM 

Removed  commented-out  code. 

Removed  unused  FORMAT  statements. 

Updated  the  version  number  and  version  date. 

RTA.FOR 

Subroutine  RTA 

Changed  argument  MLT  to  MLON. 

Variables  MLAT  and  MLON  are  now  passed  to  routine  CORRECT  1  instead  of  variables  LATADJ 
and  LTADJ. 

Removed  local  variables  LATADJ  and  LTADJ  since  they  are  no  longer  used. 

Removed  commented-out  code. 

Subroutine  COR_MAX 

Replaced  variable  EE  (molecular  ion  density  at  height  hmE)  with  argument  EE+FE  (total  ion  density 
at  height  hmE)  in  the  first  call  to  routine  GET_ONEJFO  so  that  the  correct  ratio  nmEnew/nmEoid 
(variable  FRE)  is  returned. 

Removed  variables  NME  and  DNE  from  the  first  call  to  routine  GET_ONE__FO. 

Replaced  variable  FF  (0+  density  at  height  hmF2)  with  argument  EF+FF  (total  ion  density  at  height 
hmF2)  in  the  second  call  to  routine  GET_ONEJFO  so  that  the  correct  ratio  nmF2,new/nmF2,0id  (variable 
FRF)  is  returned. 

Removed  variables  NMF2  and  DNF  from  the  second  call  to  routine  GET_ONE_FO. 

Variables  EE,  EF,  FE,  and  FF  are  now  passed  to  routine  DO__ADJ  instead  of  variables  DNE,  DNF, 
and  DNM. 

Removed  local  variables  NME,  NMF2,  DNE,  DNF,  and  DNM  since  they  are  no  longer  used. 

Removed  commented-out  code. 

Subroutine  DO_ADJ 

Complete  rewrite  to  change  algorithm  and  to  improve  coding  and  internal  documentation. 

E-layer  and  F2-layer  density  adjustments  are  now  done  on  both  the  0+  and  molecular  ion  density 
profiles. 

Subroutine  GET_ONE_FO 

Removed  variables  NMAX  and  DNM  AX  from  argument  list. 

Removed  variables  NMAX  and  DNMAX  since  they  are  no  longer  used. 

Removed  commented-out  code. 

Removed  unused  FORMAT  statements. 

- 

TECCALC.FOR 

Subroutine  LOCATE 

Modified  logic  to  eliminate  .EQV.  logical  operator,  which  caused  a  run-time  access  violation  error 
under  Microsoft  Fortran  PowerStation. 
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MEMORANDUM 


DATE: 

3  O-September-1996 

TO: 

Rob  Daniell 

FROM: 

Lincoln  Brown 

RE: 

Changes  to  PRISM  1. 6b  for  PRISM  1. 7 

The  changes  to  PRISM  1.6b  for  PRISM  1.7  focus  on  replacing  the  LLF  coefficient  set  again,  on  improving  the 

handling  of  collocated  data  in  the  real-time  adjustment,  and  on  improving  the  merging  of  the  parameterized  model 

density  profiles.  The  changes  are  summarized  as  follows: 

I.  BUG  FIX:  The  LLF  parameterized  model  coefficient  set  has  been  replaced  again.  Because  of  a  bug  in  the 
processing  of  the  LOWLAT  output,  the  magnetic  latitude  grid  in  the  PRISM  1.6b  LLF  coefficient  files  was 
defined  to  start  at  -33°N  instead  of  the  correct  value  of -34°N.  This  resulted  in  an  erroneous  northward  1° 
magnetic  latitude  shift  in  the  0+  densities  from  the  LLF  parameterized  model  The  LLF  coefficient  files  have 
been  regenerated  with  the  correct  starting  value  for  the  magnetic  latitude  grid.  Note  that  this  fix  requires  no 
change  to  the  PRISM  source  code. 

II.  BUG  FIX:  The  conversion  of  UT  from  hours  to  HHMM  in  routine  PARAM  for  the  URSI  foF2  model  has 
been  corrected.  Previously,  the  conversion  of  UT  from  decimal  hours  to  HHMM  resulted  in  truncation  of  the 
minutes  to  zero,  e.g.  a  UT  of  1.5  hours  was  erroneously  converted  to  0100.  The  conversion  has  been  modified  so 
that  the  minutes  are  not  truncated  to  zero,  e.g.  a  UT  of  1.5  hours  is  correctly  converted  to  0130. 

HI.  Based  on  work  done  by  Pat  Doherty,  we  now  recommend  that  a  27-day  running  mean  of  F  10.7  be  used 

instead  of  a  daily  value.  Pat  found  that  a  mean  Fi0.7  correlates  much  better  with  the  ionosphere  than  a  daily  value, 
either  because  the  solar  EUV  has  less  day-to-day  variability  than  F10.7  or  because  the  day-to-day  variability  in 
solar  EUV  is  uncorrelated  with  the  day-to-day  variability  in  F10.7*  The  same  applies  to  Sunspot  Number  (use  a 
27-day  running  mean  instead  of  a  daily  value),  although  we  recommend  using  F10.7  instead  of  Sunspot  Number  as 
the  solar  activity  index. 

IV.  Multiple  data  points  critically  close  to  an  output  grid  point  are  now  handled  correctly  in  the  real-time 
adjustment.  Previously,  if  one  or  more  data  points  was  within  the  critical  distance  of  an  output  grid  point,  then 
the  correction  returned  by  the  real-time  adjustment  was  the  correction  from  the  last  data  point  within  the  critical 
distance.  Now,  if  one  or  more  data  points  are  within  the  critical  distance  of  an  output  grid  point,  then  the 
correction  returned  from  the  real-time  adjustment  is  the  average  of  corrections  from  those  data  points. 

V.  The  top-level  parameterized  model  routine  PARAM  has  undergone  substantial  changes.  The  original 
motivation  was  a  report  from  Jim  Secan  at  Northwest  Research  Associates  illustrating  a  discontinuity  in  PRISM 
1.6b  TEC  at  the  latitude  transition  region  between  the  mid-latitude  and  high-latitude  parameterized  models.  A 
closer  examination  of  routine  PARAM  has  yielded  the  following  improvements: 

A.  It  has  been  completely  rewritten  for  optimization  and  readability.  The  number  of  routine  calls  has  been 
reduced. 

B.  The  transition  between  the  mid-latitude  and  high-latitude  parameterized  models  is  smoother. 

Previously,  in  the  mid/high-latitude  transition  region,  the  mid-latitude  peak  density  and  height  were  used  to 
merge  the  mid-latitude  and  high-latitude  density  profiles,  resulting  in  a  discontinuity  in  density  at  the 
boundary  between  the  high-latitude  region  and  the  mid/high-latitude  transition  region.  Now,  weighted 
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averages  of  mid-latitude  and  high-latitude  peak  densities  and  heights  are  used  in  the  merging,  resulting  in  a 
smoother  transition. 


C.  The  transition  between  the  low-latitude  and  mid-latitude  parameterized  models  is  smoother.  Previously, 
a  weighted  average  of  low-latitude  and  mid-latitude  critical  frequencies  was  used  to  calculate  the  peak  density 
for  the  merging  of  the  low-latitude  and  mid-latitude  density  profiles.  Now,  the  peak  density  used  in  the 
merging  is  calculated  from  a  weighted  average  of  low-latitude  and  mid-latitude  peak  densities.  This  is 
consistent  with  the  profile  merging  process,  and  results  in  a  smoother  transition  between  the  low-latitude  and 
mid-latitude  parameterized  models. 

VI.  The  mid-level  parameterized  model  routines  LOW__PARAM,  MIDJPARAM,  and  USUMODEL  have  been 
rewritten  to  accommodate  changes  in  the  top-level  parameterized  model  routine  PARAM,  to  provide  small 
performance  gains  by  reducing  routine  calls  and  eliminating  unnecessary  calculations,  and  to  improve  readability. 

VII.  The  low-level  parameterized  USU  model  routines  RECON,  FOUR_COEFF,  and  FULLJDATA  have  been 
modified  to  accommodate  changes  to  the  mid-level  parameterized  model  routine  USUMODEL. 

Vin.  Several  routines  in  module  ENVIRON.FOR  have  been  replaced  for  compatibility  with  PIM. 

IX.  Several  routines  have  been  modified  to  accommodate  the  new  version  of  routine  SOLDEC. 

X.  Several  minor  changes  have  been  made  for  compatibility  with  Microsoft  FORTRAN.  They  do  not  impact 
the  results. 

The  table  below  describes  the  changes  that  I  made  to  PRISM  1.6b  to  produce  PRISM  1.7. 
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Module 

Program  Unit 

Description  of  Changes 

ADJ._FNH.FOR 

Subroutine  GET__FOS 

Removed  argument  ERR  from  the  call  to  INFLECTION  since  it  is  no  longer  used  by  that  routine. 
Removed  local  variable  ERR  since  it  is  no  longer  used. 

Subroutine  INFLECTION 

Removed  input  argument  ERR  since  it  is  not  used. 

Changed  declaration  of  local  variable  TEST  from  LOGICAL*  1  to  LOGICAL. 

Subroutine  GET_FOE 

Removed  since  it  is  no  longer  used. 

ENVIRON.  FOR 

Subroutine  SOLZA 

Removed  since  it  has  been  replaced  by  routine  SOLSZA. 

Subroutine  SSOLPT 

Removed  since  it  has  been  replaced  by  routine  SOLSUB. 

Function  SOLDEC 

Replaced  by  routine  of  the  same  name. 

Function  SOLANG 

New  routine. 

Subroutine  SOLSUB 

New  routine. 

Function  SOLSZA 

New  routine. 

FMODEL.FOR 

Subroutine  FMODEL 

Removed  argument  HBOT  from  call  to  PARAM  since  it  is  not  used  by  that  routine. 

Removed  last  argument  in  call  to  PARAM  since  it  is  no  longer  used  by  that  routine. 

Removed  local  variable  HBOT  since  it  is  no  longer  used. 

GETDAT.FOR 

Subroutine  GETDAT 

Removed  argument  ESWITCH  from  call  to  MODEL  FLAGS  since  it  is  no  longer  used. 

Subroutine  MODELJFLAGS 

Removed  assignment  of  input  argument  ESW. 

Removed  input  argument  ESW  since  it  is  no  longer  used. 

Subroutine  PREC__DATA 

Changed  declaration  of  local  variable  TEST  from  LOGICAL*  1  to  LOGICAL. 

Subroutine  J4ARR 

Changed  declaration  of  local  variable  TEST  from  LOGICAL*  1  to  LOGICAL. 

HLIM.FOR 

Subroutine  REGMOD 

Removed  argument  HBOT  from  call  to  PARAM  since  it  is  not  used  by  that  routine. 

Removed  last  argument  in  call  to  PARAM  since  it  is  no  longer  used  by  that  routine. 

Removed  local  variable  HBOT  since  it  is  no  longer  used. 

Replaced  call  to  SOLZA  with  call  to  SOLSZA. 

HLISM.FOR 

Subroutine  HLISM 

Removed  local  variable  MFLAG  since  it  is  no  longer  used. 

ENDIRECT.INC 

n/a 

Removed  common  block  INDIRECT  variable  ES WTTCH  since  it  is  no  longer  used. 

INIT.FOR 

Subroutine  LDITER 

Changed  declaration  of  output  argument  LAYFLG  from  LOGICAL*  1  to  LOGICAL. 

Subroutine  GETFC 

Changed  declaration  of  local  variable  ERROR  from  L0GICAL*1  to  LOGICAL. 

!  Subroutine  LUDCMP 

Changed  PAUSE  statement  to  STOP  statement. 

IOJJTIL.FOR 

Subroutine  UDETID 

j 

Added  input  argument  UT. 

Argument  DAY  is  no  longer  converted  to  a  real  number  in  call  to  SOLDEC. 

Added  argument  UT  to  call  to  SOLDEC. 

LAYFLG.INC 

n/a 

Changed  declaration  of  common  block  LAYFLG  variable  LAYFLG  from  LOGICAL*  1  to 

LOGICAL. 
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Module 


LOW_P  ARA.  FOR  Subroutine  LOW_PARAM 


MATH_UTI.FOR  Subroutine  1NTRP 

Subroutine  SVDCMP 


MID_PARA.FOR  Subroutine  MID_P  ARAM 


OUTPUT.FOR  I  Subroutine  FINAL 


PARAM.FOR  Subroutine  PARAM 


PRISM.FOR 


RE  ADJDB  A.  FOR 


USERJNP.FOR 


Program  PRISM 


Subroutine  READ  DBASES 


Subroutine  RDLOW 


Subroutine  LDETID 


Subroutine  READ_E 


Subroutine  EDETID 


Subroutine  READMID 


Subroutine  DETID 


Subroutine  READUSU 


Subroutine  RTA 


Subroutine  CORRECT  1 


Subroutine  COR  MAX 


Subroutine  GET  ONE  FO 


Subroutine  COR  PRO 


Subroutine  CHECKSTAT 


Subroutine  GIVE  DATA 


Description  of  Changes  ( continue 


Rewrite  for  optimization  and  readability. 

Removed  input  argument  LAYR  since  it  is  no  longer  used. 
Changed  intrinsic  function  ALOG  to  its  generic  equivalent  LOG. 


Modified  logic  to  eliminate  arithmetic  IF  statements. 


Changed  PAUSE  statement  to  PRINT  statement. 


Rewrite  for  optimization  and  readability. 

Removed  input  argument  LAYR  since  it  is  no  longer  used. 
Changed  intrinsic  function  ALOG  to  its  generic  equivalent  LOG. 


Changed  declaration  of  local  variable  TEST  from  LOGICAL*!  to  LOGICAL. 


Rewrite  for  optimization  and  readability. 

Fixed  a  bug  in  the  conversion  of  UT  from  hours  to  HHMM  for  the  call  to  F2URSI. 

In  the  mid/high-latitude  transition  region,  weighted  averages  of  mid-latitude  and  high-latitude  peak 
densities  and  heights  are  now  used  in  the  merging  of  profiles  instead  of  the  mid-latitude  peak  density 
and  height. 

In  the  low/mid-latitude  transition  region,  the  peak  density  used  to  merge  profiles  is  calculated  from  a 
weighted  average  of  low-latitude  and  mid-latitude  peak  density  instead  of  a  weighted  average  of  low- 
latitude  and  mid-latitude  critical  frequency,  to  be  consistent  with  the  profile  merging  process. 
Simplified  logic  since  LAYR  flag  is  no  longer  used. 

Removed  input  argument  LAYR  since  it  is  no  longer  used. 

Removed  input  argument  HBOT  since  it  is  not  used. 

Removed  last  argument  from  calls  to  USUMODEL  since  it  is  no  longer  used. 

A  single  call  to  USUMODEL  is  now  used  instead  of  two  calls. 

Removed  last  argument  from  calls  to  MID_PARAM  since  it  is  no  longer  used. 

A  single  call  to  MID_PARAM  is  now  used  instead  of  two  calls. 

Removed  last  argument  from  calls  to  LOW_PARAM  since  it  is  no  longer  used. 

A  single  call  to  LOW_PARAM  is  now  used  instead  of  two  calls. 

Removed  call  to  GET_FOE  since  it  is  no  longer  used. 

Removed  logic  involving  common  block  INDIRECT  variable  ESWITCH  since  it  is  no  longer  used. 
Changed  intrinsic  function  ALOG  to  its  generic  equivalent  LOG. 


Undated  the  version  number  and  version  date. 


Added  argument  UT  to  call  to  RDLOW. 


Added  input  argument  UT. 

Added  argument  UT  to  call  to  LDETID. 


Added  input  argument  UT. 

Argument  DAY  is  no  longer  converted  to  a  real  number  in  call  to  SOLDEC. 
Added  argument  UT  to  call  to  SOLDEC. 


Added  argument  UT  to  call  to  EDETID. 


Added  input  argument  UT. 

Argument  DAY  is  no  longer  converted  to  a  real  number  in  call  to  SOLDEC. 
Added  argument  UT  to  call  to  SOLDEC. 


Added  argument  UT  to  call  to  DETID. 


Added  input  argument  UT. 

Argument  DAY  is  no  longer  converted  to  a  real  number  in  call  to  SOLDEC. 
Added  argument  UT  to  call  to  SOLDEC. 


Added  argument  UT  to  calls  to  UDETID. 


Changed  declaration  of  local  variable  IERR  from  LOGICAL*!  to  LOGICAL. 


Rewrite  for  optimization  and  readability. 

If  one  or  more  data  points  is  within  the  critical  distance  of  the  point  of  interest,  then  the  correction  at 
the  point  of  interest  is  now  calculated  as  the  average  of  corrections  from  those  data  points  instead  of 
the  correction  of  the  last  data  point  within  the  critical  distance  being  used. 


Changed  declaration  of  output  argument  IERR  from  LOGICAL*!  to  LOGICAL. 


Changed  declaration  of  output  argument  IERR  from  LOGICAL*!  to  LOGICAL. 


Changed  declaration  of  output  argument  IERR  from  LOGICAL*!  to  LOGICAL. 


Changed  declaration  of  input  argument  TEST  from  LOGICAL*!  to  LOGICAL. 


Changed  declaration  of  local  variable  TEST  from  LOGICAL*!  to  LOGICAL. 
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Module 

Program  Unit 

Description  of  Changes  {continued) 

USUMODEL.FOR 

Subroutine  USUMODEL 

Rewrite  for  optimization  and  readability. 

Removed  input  argument  LAYR  since  it  is  no  longer  used. 

Changed  intrinsic  function  ALOG  to  its  generic  equivalent  LOG. 

Subroutine  RECON 

Replaced  input  argument  XNMLAT  with  input  arguments  SMLATE  and  SMLATF. 

Replaced  argument  XNMLAT  with  argument  SMLATE  and  SMLATF  in  call  to  FOUR_COEFF. 
Removed  argument  LAYR  from  calls  to  FOUR_COEFF  and  FULLJDATA. 

Removed  input  argument  LAYR  since  it  is  no  longer  used.  ■ 

Subroutine  FOUR__COEFF 

Replaced  input  argument  XNMLAT  with  input  arguments  SMLATE  and  SMLATF. 

Fourier  coefficients  are  now  calculated  for  all  three  ions. 

Removed  input  argument  LAYR  since  it  is  no  longer  used. 

Subroutine  FULL.DATA 

Altitude  profiles  are  now  calculated  for  all  three  ions. 

Removed  input  argument  LAYR  since  it  is  no  longer  used. 

The  figure  below  illustrates  the  improvement  in  the  merging  of  the  regional  parameterized  models.  It  shows  a  profile 
of  foF2  vs.  geographic  latitude.  Note  the  discontinuity  in  PRISM  1.6b  f<>F2  in  the  mid/high-latitude  transition  region 
and  the  smoother  transition  in  the  PRISM  1.7  profile.  The  differences  at  low  latitudes  are  due  to  the  change  in  the 
low/mid-latitude  transition  method  and  the  corrected  LLF  parameterized  model  coefficients. 
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test. out 

Year:  1996 
Day:  80 
UT:  12.00  hr 

F10.7:  193.0 
SSN:  150.0 
Kp:  3.0 

Lon:  150.0  °E  geographic 
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MEMORANDUM 


DATE: 

7 -November- 199  6 

TO: 

Rob  Daniell 

FROM: 

Lincoln  Brown 

RE: 

Changes  to  PRISM  1. 7  for  PRISM  1. 7a 

The  changes  to  PRISM  1.7  for  PRISM  1.7a  focus  on  improvements  to  the  real-time  data  ingestion.  The  changes  are 
summarized  as  follows: 

I.  BUG  FIX:  The  presence  of  all  five  possible  DISS  real-time  data  files  no  longer  causes  an  array-out-of- 
bounds  error  in  the  DISS  data  ingestion.  Previously,  if  root  names  for  all  five  DISS  data  files  were  given  in 
the  PATH__NAM.TXT  file,  then  the  DISS  data  ingestion  would  attempt  to  access  a  nonexistent  sixth  element 
of  the  array  that  holds  the  DISS  data  file  root  names.  This  problem  has  been  eliminated  due  to  the  rewrite  of 
the  DISS  data  ingestion  (see  item  5). 

II.  BUG  FIX:  The  presence  of  all  five  possible  IMS  real-time  data  files  no  longer  causes  an  array-out-of- 
bounds  error  in  the  IMS  data  ingestion.  Previously,  if  root  names  for  all  five  IMS  data  files  were  given  in  the 
PATH_NAM.TXT  file,  then  the  IMS  data  ingestion  would  attempt  to  access  a  nonexistent  sixth  element  of 
the  array  that  holds  the  IMS  data  file  root  names.  This  problem  has  been  eliminated  due  to  the  rewrite  of  the 
IMS  data  ingestion  (see  item  6). 

HI.  BUG  FIX:  The  presence  of  all  eight  possible  DMSP  real-time  data  files  no  longer  causes  an  array-out- 
of-bounds  error  in  the  DMSP  data  ingestion.  Previously,  if  root  names  for  all  eight  DMSP  data  files  were 
given  in  the  PATH_NAM.TXT  file,  then  the  DMSP  data  ingestion  would  attempt  to  access  a  nonexistent 
sixth  element  of  the  array  that  holds  the  DMSP  data  file  root  names.  This  problem  has  been  eliminated  due 
to  the  rewrite  of  the  DMSP  data  ingestion  (see  item  7). 

IV.  BUG  FIX:  Large  amounts  of  DISS  foF2  real-time  data  and/or  large  amounts  of  IMS  real-time  data  no 
longer  cause  an  array-out-of-bounds  error  in  the  midlatitude  correction  algorithm.  Previously,  array 
allocation  in  the  midlatitude  correction  algorithm  for  corrections  for  these  kinds  of  data  was  inconsistent  with 
limits  defined  in  the  real-time  data  ingestion.  This  problem  was  reported  by  David  Coxwell  at  Air  Force 
Institute  of  Technology.  This  problem  has  been  eliminated  due  to  the  rewrite  of  the  DISS  and  IMS  data 
ingestion  (see  items  5  and  6)  and  the  rewrite  of  the  midlatitude  correction  algorithm  (see  item  9). 

V.  The  DISS  real-time  data  ingestion  has  been  rewritten  for  speed  and  readability.  A  single  pass  through 
the  DISS  data  is  now  required  instead  of  two  passes.  Up  to  850  DISS  data  records  can  now  be  accepted  (i.e. 
meet  criteria  for  use  in  the  real-time  adjustment)  from  all  DISS  data  files  combined,  a  reduction  from  the 
previous  limit  of  7000  accepted  records  per  DISS  data  file.  The  limit  of  850  accepted  DISS  data  records  is 
based  on  50  DISS  stations  and  a  15  minute  time  resolution  in  a  [-2,2]  hour  time  window  relative  to  the  time 
of  the  run. 

VI.  The  IMS  data  ingestion  has  been  rewritten  for  speed  and  readability.  A  single  pass  through  the  IMS 
data  is  now  required  instead  of  two  passes.  Up  to  21600  IMS  data  records  can  now  be  accepted  (i.e.  meet 
criteria  for  use  in  the  real-time  adjustment)  from  all  IMS  data  files  combined,  a  reduction  from  the  previous 
limit  of  7000  accepted  records  per  IMS  data  file.  The  limit  of 21600  accepted  IMS  data  records  is  based  on 
200  IMS  stations  reporting  12  simultaneous  measurements  and  a  time  resolution  of  15  minutes  in  a  [-2,2] 
hour  time  window  relative  to  the  time  of  the  run. 
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VII.  The  DMSP  data  ingestion  has  been  rewritten  for  speed  and  readability.  A  single  pass  through  the 
DMSP  data  is  now  required  instead  of  three  passes.  Up  to  24  sets  of  each  of  the  three  kinds  of  DMSP  data 
are  now  accepted  by  PRISM  from  all  DMSP  data  files  combined,  a  reduction  from  the  previous  limit  of  24 
sets  per  DMSP  data  file.  The  limit  of  24  sets  is  based  on  8  satellites,  an  orbital  period  of  6060  seconds,  and 
a  [-2,2]  hour  time  window  relative  to  the  time  of  the  run.  Up  to  1451  SSIES  ION  DRIFT  data  records  can  be 
accepted  (i.e.  meet  criteria  for  use  in  determining  equatorward  trough  boundaries)  per  set,  a  limit  which  has 
not  changed.  The  limit  of  145 1  accepted  SSIES  ION  DRIFT  data  records  is  based  on  a  5  second  time 
resolution  in  a  [-2,0]  hour  time  window  relative  to  the  time  of  the  run,  with  10  extra  records  allowed  for 
orbital  overlap.  Up  to  3848  SSIES  IN  SITU  PLASMA  data  records  can  now  be  accepted  (i.e.  meet  criteria 
for  use  in  the  real-time  adjustment)  from  all  DMSP  data  files  combined,  a  reduction  from  the  previous  limit 
of  145 1  accepted  records  per  DMSP  data  file.  The  limit  of 3848  accepted  SSIES  IN  SITU  PLASMA  data 
records  is  based  on  8  satellites  and  a  time  resolution  of  30  seconds  in  a  [-2,2]  hour  time  window  relative  to 
the  time  of  the  run.  Up  to  721 1  SSJ/4  data  records  can  be  accepted  (i.e.  meet  criteria  for  use  in  determining 
auroral  oval  boundaries  and  for  use  in  the  real-time  adjustment)  per  set,  a  limit  which  has  not  changed.  The 
limit  of  721 1  accepted  SSJ/4  data  records  is  based  on  a  time  resolution  of  1  second  in  a  [-2,0]  hour  time 
windows  relative  to  the  time  of  the  run,  with  10  extra  records  allowed  for  orbital  overlap. 

VHI.  The  output  station  list  determination  has  been  rewritten  for  speed  and  readability.  A  single  pass 
through  the  output  station  list  is  now  required  instead  of  two  passes.  The  number  of  output  stations  is  now 
unlimited,  a  change  from  the  previous  limit  of  7000  stations. 

IX.  The  midlatitude  correction  algorithm  has  been  rewritten  for  the  following  reasons: 

A.  To  optimize  for  speed. 

B.  To  prevent  extraneous  zero  corrections  from  being  included  in  the  real-time  adjustment.  The 
extraneous  corrections  were  not  based  on  valid  data,  and,  although  zero,  influenced  the  correction  fields 
and  increased  run  time. 

C.  To  remove  coding  for  BOTTOMSIDE  data  since  that  data  type  is  no  longer  considered  by  PRISM. 

D.  To  remove  redundant  and  commented-out  coding. 

E.  To  improve  readability. 

X.  Limits  on  the  ingestion  of  real-time  data  are  now  globally  defined  in  new  INCLUDE  file  rtdlimit.inc. 

XI.  The  HLE  model  has  been  removed  since  it  is  no  longer  used.  All  components  of  PRISM  upon  which 
HLE  solely  depends  have  also  been  removed,  such  as  the  MSIS-86  neutral  atmosphere  model  and  a  number 
of  INCLUDE  files.  The  HLE  data  files  CHEM.FIL,  TIMING.FIL,  and  SOLAR.FIL  are  no  longer  used.  The 
path  to  the  HLE  files  is  still  read  from  the  PATH_NAM.TXT  file  for  compatibility,  but  it  is  not  used. 

XII.  Coding  for  the  BOTTOMSIDE  data  type  has  been  removed  since  it  is  no  longer  considered  in  PRISM. 
Xni.  INCLUDE  file  dlta.inc  has  been  removed  since  it  is  not  used. 

XIV.  A  number  of  unused  PARAMETERS  and  variables  have  been  removed. 

The  performance  gain  due  to  these  changes  is  hard  to  predict  since  it  depends  on  the  amount  and  character  of  the 
real-time  data.  For  many  of  the  standard  test  cases  provided  with  delivery,  run-time  has  been  decreased  by  about 
15%.  The  memory  requirement  has  been  reduced  by  about  0.5  MB,  a  7%  decrease. 

The  table  below  describes  the  changes  that  I  made  to  PRISM  1.7  to  produce  PRISM  1.7a. 
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Module 


ALPHAN.INC 


CHEM.INC 


DLTA.INC 


ESPECT.INC 


CGM  UTIL.FOR 


ENVIRON.  FOR 


GETDAT.FOR 


Program  Unit 


n/a 


n/a 


n/a 


n/a 


Subroutine  BNDRY 


Subroutine  NEUATM 


Subroutine  GTS5 


Function  DENSS 


Function  GLOBE5 


Function  GLOB5L 


Function  DNET 


Function  CCOR 


Block  Data  PRMDTD 


Function  SOLANG 


Subroutine  SOLSUB 


Function  SOLSZA 


Subroutine  GETDAT 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


DO_DIR 


DO JONO 


DO  TEC 


DOJDMSP 


INRTD 


INDISS 


INIMS 


INDMSP 


DETOSL 


PREC  DATA 


Subroutine  LOADDIR 


HLEMODEL.FOR  I  n/a 


HLIM.FOR  Subroutine  REGMOD 


HLISM.FOR  I  Subroutine  HLISM 


Subroutine  MATRIX 


Subroutine  INITPR 


Description  of  Changes 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Added  INCLUDE  statement  for  INCLUDE  file  rtdlimit.inc. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Replaced  call  to  routine  DO_DIR  with  calls  to  routine  INRTD  and  routine  DETOSL. 
Removed  output  parameter  SUSI  since  it  is  not  used. 

Removed  PARAMETERS  NSTA  and  MAXTYPE  since  they  are  no  longer  used. 
Removed  local  variable  NDAT  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used  (it  has  been  replaced  by  routine  INRTD 


Removed  since  it  is  no  longer  used  (it  has  been  replaced  by  routines  INDISS  and  DETOSL 


Removed  since  it  is  no  longer  used  (it  has  been  replaced  by  routine  INIMS 


Removed  since  it  is  no  longer  used  (it  has  been  replaced  by  routine  INDMSP 


New  routine. 


New  routine. 


New  routine. 


New  routine. 


New  routine. 


Output  parameter  NREC  is  now  initialized  to  zero. 

PRECIPITATION  records  are  now  written  to  data  file  at  the  current  file  position  instead  of  the  end  of 
the  file. 

A  PRECIPITATION  data  header  is  no  longer  written  to  the  data  file. 

Standard  deviation  placeholders  are  no  longer  included  in  the  data  record  written  to  the  data  file. 
INCLUDE  statements  for  INCLUDE  files  dpath.inc  and  direct_d.inc  have  been  removed  since  they 
are  no  longer  needed. 

Removed  local  variable  DUMMY  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used. 


Removed  coding  for  HLE  model  since  it  is  no  longer  used. 

Removed  input  parameters  F10ANA,  RF10NA,  and  APNA  since  they  are  no  longer  used. 
Removed  local  variables  SZA,  HMED,  and  ZNEU  since  they  are  no  longer  used. 

Removed  PARAMETER  NZNEU  since  it  is  no  longer  used. 

Removed  coding  for  BOTTOMSIDE  data  since  the  BOTTOMSEDE  data  type  is  no  longer  used. 
Removed  INCLUDE  statement  for  INCLUDE  file  phys4  co.inc  since  it  is  not  needed. 


Removed  arguments  RF10P7,  F10P7A,  and  AP  from  call  to  routine  MATRIX  since  they  are  no 
longer  used  by  that  routine. 

Removed  input  parameters  RF10P7,  F10P7A,  and  AP  since  they  are  no  longer  used. 


Removed  arguments  F10P7A,  RF10P7,  and  AP  from  calls  to  routine  REGMOD  since  they  are  no 
longer  used  by  that  routine. 

Removed  input  parameters  RFI0P7,  F10P7A,  and  AP  since  they  are  no  longer  used. 


Added  INCLUDE  statement  for  INCLUDE  file  rtdlimit.inc. 
PARAMETER  MAXORB  has  been  renamed  to  MORBIT. 


Removed  since  it  is  no  longer  used. 


Removed  PARAMETERS  LUIDL,  LUNATM,  LUNUMB,  LURATE,  LURI,  LURO,  LUSGP,  and 
LUTEMP  since  they  are  not  used. 

Removed  PARAMETERS  LUCHEM,  LUSOLR,  and  LUTIME  since  they  are  no  longer  used. _ 


Removed  since  it  is  no  longer  used. 
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Module 


MATHJJTI.FOR 


MATH4  CO.INC 


MIDLAT.FOR 


Program  Unit 


Subroutine  RTBIS 


Subroutine  INTRP 


Subroutine  CHMRRC 


n/a 


Subroutine  MIDLAT 


NATMOS.INC 


NEUATM..INC 


NEWFIT.INC 


Subroutine  PER_j\RR 


n/a 


n/a 


n/a 


OUTPUT.FOR  I  Subroutine  GRIDJDUTPUT 


Subroutine  STATION„OUTPUT 


PHANTOM. FOR  1  Subroutine  PHANTM 


PHYS4_CO.INC 


PRECIP.FOR 


PRECIP.INC 


PRISM.FOR 


n/a 


Subroutine  PRECIP 


Block  Data  PREC 


Program  PRISM 


PSPECT.INC  I  _ n/a _ 


READ JDBA.FOR  Subroutine  DBASES 

Subroutine  READ__DBASES 

Subroutine  RDFILES _ 

Subroutine  RPC  HEM _ 

Subroutine  RDTIME _ 

Subroutine  RDSOLR 


Description  of  Changes  (continue 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used.  _ 


Removed  PARAMETERS  DPH,  PIQ2,  and  TWOPI  since  they  are  not  used. _ 


Rewrite  for  the  following  reasons: 

1.  To  optimize  calculation  of  corrections. 

2.  To  prevent  extraneous  zero  corrections  from  being  included  in  the  real-time  adjustment.  The 
extraneous  corrections  were  not  based  on  valid  data,  and  although  zero,  did  influence  the 
correction  fields,  and  increased  run  time. 

3.  To  remove  coding  for  BOTTOMSIDE  data  since  that  data  type  is  no  longer  used. 

4.  To  remove  redundant  and  commented-out  coding. 

5.  To  improve  readability. 

Removed  arguments  F10P7A,  RF10P7,  and  AP  from  calls  to  routine  REGMOD  since  they  are  no 
longer  used  by  that  routine.  _ 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used.  _ 


Combined  common  blocks  INWFT  and  NEWFT  into  common  block  NEWFIT. 

Removed  variables  NHFM,  NFFM,  NHEM,  NFEM,  NB1M,  NB1D,  NB2M,  NB2D,  NHTM,  NNTM, 
FFDATA,  EFDATA,  FHDATA,  EHDATA,  BIDATA,  B2DATA,  NTDATA,  HTDATA,  BT1LAT, 
BT1LON,  BT2LAT,  BT2LON,  B1C,  and  B2C  from  common  block  NEWFIT  since  they  are  no 
longer  used. 

Replaced  PARAMETERS  NM AX  and  NNMAX  with  PARAMETERS  MFFD,  MHFD,  MFED, 
MHED,  MNTD,  and  MHTD.  The  new  PARAMETERS  are  dependent  on  PARAMETERS  in 
INCLUDE  file  rtdlimit.inc.  _ 


Removed  arguments  F10P7A,  RF10P7,  and  AP  from  call  to  routine  REGMOD  since  they  are  no 
longer  used  by  that  routine. 

Removed  INCLUDE  statement  for  INCLUDE  file  dlta.inc  since  it  is  not  used. _ 


Removed  arguments  F10P7A,  RF10P7,  and  AP  from  call  to  routine  REGMOD  since  they  are  no 
longer  used  by  that  routine. 


Removed  arguments  F10P7A,  RF10P7,  and  AP  from  calls  to  routine  REGMOD  since  they  are  no 
longer  used  by  that  routine. 

Removed  INCLUDE  statement  for  INCLUDE  file  indirect.inc  since  it  is  no  longer  needed. 


Removed  since  it  is  no  longer  used.  _ 


Added  INCLUDE  statement  for  INCLUDE  file  rtdlimit.inc. _ 


PARAMETER  MAXORB  has  been  moved  to  INCLUDE  file  rtdlimit.inc  and  renamed  to  MORBIT. 


Added  INCLUDE  statement  for  INCLUDE  file  rtdlimit.inc. 

PARAMETER  MAXORB  has  been  renamed  to  MORBIT. _ 


Removed  argument  SUSI  from  call  to  routine  GETDAT  since  it  is  not  used. 

Removed  arguments  RF10P7,  F10P7A,  and  AP  from  calls  to  routine  HLISM  since  they  are  no  longer 
used  by  that  routine. 

Removed  INCLUDE  statement  for  INCLUDE  file  dlta.inc  since  it  is  not  used. 

Undated  the  version  number  and  version  date.  . . . . . 


Removed  since  it  is  no  longer  used.  _ 


Removed  argument  PLHE(J9:J10)  from  call  to  routine  READJDBASES  since  it  is  no  longer  used  by 
that  routine. 

Removed  local  variables  J9  and  JO  since  they  are  no  longer  used. _ 


Removed  call  to  routine  RDFILES  since  it  is  no  longer  needed. 

Removed  input  parameter  HLEPATH  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used.  _ 


Removed  since  it  is  no  longer  used. 


Removed  since  it  is  no  longer  used.  ____ 


Module 

Program  Unit 

Description  of  Changes  ( continued) 

RTA.FOR 

Subroutine  RTA 

Added  INCLUDE  statement  for  INCLUDE  file  rtdlimit.inc. 

Removed  arguments  FFDATA,  EFDATA,  FHDATA,  EHDATA,  NTDATA,  and  HTDATA  from 
calls  to  routine  CORRECT1  since  they  are  no  longer  used  by  that  routine. 

Removed  arguments  NMAX  and  NNMAX  from  calls  to  routine  CORRECT  1  since  they  are  no 
longer  used  by  that  routine. 

Removed  coding  for  BOTTOMSIDE  data  since  the  BOTTOMSIDE  data  type  is  no  longer  used. 
Removed  INCLUDE  statement  for  INCLUDE  file  dlta.inc  since  it  is  not  used. 

Subroutine  CORRECT  1 

Validity  checking  is  no  longer  needed  since  all  data  points  reaching  this  routine  are  to  be  used. 
Removed  input  argument  MDAT  since  it  is  no  longer  used. 

Removed  input  argument  LDAT  since  it  is  no  longer  used. 

RTDLIMIT.INC 

n/a 

New  file. 

SOLAR.INC 

n/a 

Removed  since  it  is  no  longer  used. 

SPECIE.INC 

n/a 

Removed  since  it  is  no  longer  used. 

STRINGS.FOR 

Subroutine  STRTRM 

New  routine. 

TIME. INC 

n/a 

Removed  since  it  is  no  longer  used. 

MEMORANDUM 


DATE: 

14-February-1997 

m 

TO: 

Rob  Daniell 

T* 

FROM: 

Lincoln  Brown 

RE: 

Changes  to  PRISM  1. 7a  for  PRISM  1. 7b 

The  changes  to  PRISM  1 .7a  for  PRISM  1 .7b  focus  on  miscellaneous  bug  fixes.  The  changes  are  summarized  as 
follows: 

1.  BUG  FIX:  The  maximum  allowed  number  of  real-time  IMS  data  points  that  PRISM  can  accept  has  been 
changed  from  21,600  to  40,800  in  order  to  be  consistent  with  the  I/O  Specification  (“200  IMS  stations 
reporting  12  simultaneous  measurements  and  a  time  resolution  of  15  minutes  in  a  +/-  2  hour  time  window 
centered  on  the  date  and  UT  of  the  run”,  or  200  12-[l+4-60/15]=40,800).  The  previous  value  of  21,600  was 
consistent  with  a  time  resolution  of  30  minutes  (200- 1 2 •[  1  +4-60/3 0]  =2 1,600). 

2.  BUG  FIX:  Array  allocation  for  real-time  data  has  been  increased  to  allow  for  phantom  data.  Previously,  the 
presence  of  the  maximum  allowed  number  of  real-time  data  points  in  combination  with  phantom  data  points 
could  cause  an  array-out-of-bounds  error  in  the  calculation  of  the  low/midlatitude  correction  fields. 

3.  BUGFIX:  In  the  determination  of  the  output  station  list  from  accepted  real-time  DISS  data,  the  search  of 
the  direct  data  file  for  the  IONOSONDE  data  section  has  been  modified  to  prevent  a  read-past-end-of-file 
error.  In  subroutine  DETOSL  in  module  GETDAT.FOR,  variable  NSTOUT  is  now  reinitialized  to  zero 
after  a  non-IONOSONDE  data  section  in  the  direct  data  file  is  skipped  over  during  the  search.  Previously,  if 
the  direct  data  file  contained  one  or  more  non-IONOSONDE  data  sections,  but  contained  no  IONOSONDE 
data  section,  the  value  of  NSTOUT  from  the  last  non-IONOSONDE  data  section  would  be  used  to  attempt  to 
read  an  IONOSONDE  data  section  past  the  end  of  the  direct  data  file,  resulting  in  a  run-time  error. 

4.  BUG  FIX:  In  the  determination  of  the  output  station  list  from  accepted  real-time  DISS  data,  the  corrected 
geomagnetic  local  times  of  the  output  stations  now  correspond  to  the  day  of  the  year  and  Universal  Time  of 
the  run.  Previously,  the  corrected  geomagnetic  local  times  of  the  DISS  data  records  were  used  for  the  output 
stations,  resulting  in  as  much  as  a  2  hour  error  in  corrected  geomagnetic  local  time. 

PRISM’S  memory  requirement  has  been  increased  by  about  0.25  MB,  a  3%  increase.  Its  run-time  will  be  unaffected 
unless  the  output  station  list  is  taken  from  a  very  large  set  of  accepted  DISS  data. 

The  table  below  describes  the  changes  that  I  made  to  PRISM  1 .7a  to  produce  PRISM  1.7b. 
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Module 

Program  Unit 

Description  of  Changes 

GETDAT.FOR 

Subroutine  INDISS 

A  maximum  of  MDISS-MPIONO  DISS  data  records  are  now  accepted  instead  of  a  maximum  of 
MDISS  DISS  data  records  in  order  to  leave  room  for  phantom  DISS  data. 

Subroutine  INDMSP 

A  maximum  of  MISP-MPINSI SSIES  IN  SITU  PLASMA  data  records  are  now  accepted  instead  of  a 
maximum  of  MISP  SSIES  IN  SITU  PLASMA  data  records  in  order  to  leave  room  for  phantom  SSIES 
IN  SITU  PLASMA  data. 

Subroutine  DETOSL 

Variable  NOUTST  is  now  reinitialized  to  zero  after  non-IONOSONDE  data  sections  in  the  direct 
data  file  are  skipped  over. 

When  the  output  station  list  is  taken  from  the  DISS  data  list,  the  corrected  geomagnetic  local  time  of 
an  output  station  list  record  is  now  calculated  from  the  nominal  day  of  the  year  and  Universal  Time 
instead  of  being  taken  from  the  corresponding  DISS  data  record. 

PHANTOM.FOR 

Subroutine  PHANTM 

Moved  PARAMETERS  NZMLAI,  NZMLAP,  NZMLOI  and  NZMLOP  to  INCLUDE  file  rtdlimit.inc. 
Added  INCLUDE  statement  for  INCLUDE  file  rtdlimit.inc. 

PRISM.FOR 

Program  PRISM 

Updated  the  version  number  and  version  date. 

RTDLIMIT.INC 

n/a 

Moved  PARAMETERS  NZMLAI,  NZMLAP,  NZMLOI,  and  NZMLOP  from  subroutine  PHANTM. 
Added  PARAMETERS  MPIONO  and  MPINSI. 

Changed  value  of  PARAMETER  MDISS  from  850  to  850+MPIONO. 

Changed  value  of  PARAMETER  MIMS  from  21600  to  40800. 

Changed  value  of  PARAMETER  MISP  from  3848  to  3848+MPINSI. 
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M’E*M*OR*A4N‘D*U4M 


DATE:  13-January-1998 

TO:  Rob  Daniell 

FROM:  Lincoln  Brown 

RE :  Changes  to  PRISM  1. 7b  for  PRISM  1. 7c 

The  changes  to  PRISM  1.7b  for  PRISM  1.7c  focus  improving  the  midlatitude  real-time  adjustment  algorithm  and  on 

replacing  the  LLF  and  MLF  parameterized  models.  The  changes  are  summarized  as  follows: 

5.  BUGFIX:  The  calculation  of  the  topside  half-width  of  the  0+ density  has  been  corrected.  Previously,  due  to  a 
sign  convention  error,  the  topside  half-width  was  calculated  as  a  negative  value.  It  is  now  correctly  calculated  as 
a  positive  value.  Because  the  minimum  of  the  bottomside  and  topside  half-widths  is  used  by  PRISM,  and  the 
bottomside  0+  density  probably  always  falls  off  more  rapidly  than  the  topside  0+  density,  this  bug  has  probably 
never  influenced  the  results  of  PRISM. 

6.  BUG  FIX:  A  minor  bug  in  the  selection  of  magnetic  latitudes  for  interpolation  in  the  MLF  parameterized  model 
has  been  fixed.  This  bug  should  not  have  affected  past  results. 

7.  The  algorithm  for  scaling  the  molecular  ion  and  0+  model  density  profiles  using  simultaneous  E-layer  and  F-layer 
real-time  data  corrections  has  been  redesigned.  Previously,  in  the  presence  of  disparate  E-layer  and  F-layer 
corrections,  discontinuities  could  be  introduced  into  the  density  profiles,  resulting  in  unintentional  changes  in 
peak  heights  as  well  as  scaled  peak  densities  inconsistent  with  the  real-time  data  corrections.  The  new  scaling 
algorithm  preserves  the  heights  of  the  peaks  and  smoothly  transitions  the  density  profile  scaling  from  a  pure  E- 
layer  correction  at  the  molecular  ion  density  peak  height  to  a  pure  F-layer  correction  at  the  0+  density  peak 
height. 

8.  The  conversion  of  TEC  real-time  data  to  critical  frequency  corrections  is  now  based  simply  on  the  ratio  of  data 
TEC  to  model  TEC,  and  both  foF2  and  f0E  corrections  are  now  provided  from  the  conversion.  Previously,  a  more 
complicated  algorithm  based  on  the  topside  model  TEC  was  used,  and  only  foF2  corrections  were  provided. 

9.  Phantom  data  now  no  longer  influences  the  sunspot  number  iteration.  Previously,  due  to  the  placement  of  the 
call  to  routine  PHANTM  in  the  main  program,  phantom  data  was  generated  and  included  in  the  real-time  data  set 
prior  to  the  sunspot  number  iteration.  Since  the  sunspot  number  iteration  fits  an  optimum  sunspot  number  to  the 
real-time  data,  the  phantom  data  erroneously  influenced  the  outcome  of  the  sunspot  number  iteration.  The  call  to 
routine  PHANTM  has  been  moved  so  that  the  phantom  data  influences  only  the  midlatitude  real-time  adjustment. 

10.  In  order  to  eliminate  problems  in  merging  the  LLF  and  MLF  parameterized  models  due  to  differences  between 
the  two  theoretical  models  (LOWLAT  and  MIDLAT)  previously  used  as  their  basis,  the  LLF  and  MLF 
parameterized  models  have  been  regenerated  from  a  single  theoretical  model  (LOWLAT).  In  addition,  the 
following  improvements  have  been  made  to  LOWLAT: 

a.  The  equatorial  vertical  ExB  drift  and  its  radial  derivative  now  vary  smoothly  to  zero  in  the  drift  transition 
region.  Previously,  a  simple  linear  fall-off  was  used,  resulting  in  discontinuities  in  the  radial  derivative  at  the 
drift  transition  region  endpoints. 

b.  The  radial  derivative  of  the  equatorial  vertical  ExB  drift  is  now  zero  above  the  drift  transition  region. 
Previously,  it  was  never  set  to  zero  above  the  transition  region. 
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c.  The  dayside  and  nightside  electron  temperatures  now  merge  at  0600  and  1800  solar  local  time.  Previously, 
the  electron  temperature  was  discontinuous  at  0600  and  1800  solar  local  time. 

d.  The  neutral  wind  is  now  calculated  for  the  correct  geographic  longitude  at  all  points  along  a  field  line. 
Previously,  a  fixed  geographic  longitude  was  used,  resulting  in  an  error  in  the  neutral  wind  due  to  the 
magnetic  declination  of  the  field  line. 

1 1.  The  LLF  parameterized  model  has  been  extended  to  44°  in  absolute  latitude,  in  order  to  broaden  the  merge 
region  between  the  LLF  and  MLF  parameterized  models.  Previously,  the  LLF  parameterized  model  only  went  up 
to  34°  in  absolute  latitude. 

12.  The  lower  absolute  latitude  boundary  of  the  MLF  parameterized  model  is  now  34°,  in  order  to  reduce  the  error 
due  to  the  assumption  of  verticality  of  the  midlatitude  field  lines  and  to  make  sure  that  the  midlatitude  field  lines 
are  outside  the  region  of  vertical  ExB  drift.  Previously,  the  lower  absolute  latitude  boundary  of  the  MLF 
parameterized  model  was  30°. 

13.  The  merge  region  for  the  LLF  and  MLF  parameterized  models  has  been  broadened  to  the  absolute  latitude  range 
34°-44°,  in  order  to  improve  the  quality  of  the  merge.  Previously,  the  range  was  30°-34°,  too  narrow  to 
effectively  merge  the  two  parameterized  models. 

14.  Some  header  comments  have  been  corrected. 

The  plots  of  foF2  and  TEC  vs.  magnetic  latitude  below  illustrate  the  improvement  in  PRISM  1 .7c  regarding  the 
agreement  and  merging  of  the  LLF  and  MLF  parameterized  models.  Notice  the  discontinuities  in  foF2  and  TEC  at 
the  LLF/MLF  merge  region  of  PRISM  1.7b  (±30°-34°  magnetic  latitude),  and  their  absence  in  PRISM  1.7c.  The 
plots  were  generated  by  PRISM  for  the  following  conditions:  year  1981,  day  of  the  year  173,  Universal  Time  0000, 
F10.7  210,  Kp  3.5,  IMF  By  positive,  IMF  Bz  negative,  270°E  magnetic  longitude,  no  URSI  f0F2  normalization,  and  no 
real-time  data. 


The  table  below  describes  the  changes  that  I  made  to  PRISM  1.7b  to  produce  PRISM  1.7c. 
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Module 

GEN.FOR 

HLIM.FOR 

HLISM.FOR ' 
LOWER.INC 

mid_para.forT 

MIDLAT.FOR 

MIDLAT.INC 

NEWFIT.INC 

OUTPUT.FOR 

PARAM.FOR 

PH  ANTOM  .FOR 

PRISM.  FOR 

READ_DBA.FOR 

RTA.FOR 


Program  Unit 


Subroutine  GENEC 


Subroutine  REGMOD 


Subroutine  MATRIX 


Subroutine  MID_F 


Subroutine  GRID_OUTPUT 


Subroutine  STATION_OUTPUT 


Subroutine  PHANTM 


Program  PRISM 


Subroutine  LRDOPCF 


Subroutine  RTA _ 

Subroutine  COR_MAX 


Subroutine  DO__ADJ 


Subroutine  GETHW 


Description  of  Changes  


Changed  PARAMETER  MAXMX  from  40  to  50  to  accommodate  new  LLF  parameterization. _ 


Removed  calculations  of  topside  ion  column  density  and  critical  height  difference  for  TEC  data  type 
since  they  are  no  longer  needed. 

Removed  output  parameters  TCI  and  TC2  since  they  are  no  longer  used. 

Removed  calculation  of  nmF2  since  it  is  no  longer  needed. 

Removed  local  variable  NMF2  since  it  is  no  longer  needed. 


Removed  arguments  TCI  and  TC2  from  calls  to  routine  REGMOD  since  they  are  no  longer  used  by 
that  routine. 

Removed  local  variables  TCI  and  TC2  since  they  are  no  longer  used. 


Changed  PARAMETER  MOPM1  from  1 1  to  14  for  new  LLF  parameterization. 
Changed  PARAMETER  MX  from  35  to  45  for  new  LLF  parameterization. 


Corrected  test  for  magnetic  latitude  above  the  magnetic  latitude  grid  by  changing  "(AMLAT  .GT. 
ASMLAT+FLO  AT(NMLAT(I, 1  ))*  ADMLAT)"  to  "(AMLAT  .GT.  ASMLAT+FLOAT(NMLAT(I, 1  )- 
n*ADMLATV\ _ 


The  conversion  of  TEC  data  to  critical  frequency  corrections  is  now  based  solely  on  the  ratio  of  the 
data  TEC  to  the  model  TEC,  and  both  f0F2  and  f0E  corrections  are  now  provided  from  the  conversion 
process. 

Removed  arguments  TCI  and  TC2  from  calls  to  routine  REGMOD  since  they  are  no  longer  used  by 
that  routine. 

Removed  local  variables  TCI,  TC2,  W,  WFACT,  WB,  DTEC,  TFACT,  and  NMF2N  since  they  are 
no  longer  used. 

Removed  INCLUDE  statement  for  INCLUDE  file  tomid.inc  since  it  is  no  longer  needed. 


Changed  PARAMETER  MX  from  12  to  1 1  for  new  MLF  parameterization. 


Changed  value  of  PARAMETER  MFED  from  MDISS  to  MDISS+MIMS  to  accommodate  additional 
midlatitude  f0E  corrections  from  converted  TEC  data. _ . 


Removed  arguments  TCI  and  TC2  from  call  to  routine  REGMOD  since  they  are  no  longer  used  by 
that  routine. 

Removed  local  variables  TCI  and  TC2  since  they  are  no  longer  used. _ 

Removed  arguments  TCI  and  TC2  from  call  to  routine  REGMOD  since  they  are  no  longer  used  by 
that  routine. 

Removed  local  variables  TCI  and  TC2  since  they  are  no  longer  used. _ 


Changed  lower  absolute  latitude  boundary  of  pure  mid-latitude  region  from  34.  degrees  to  44. 
degrees. 

Changed  upper  absolute  latitude  boundary  of  low/mid-latitude  transition  region  from  34.  degrees  to 
44.  degrees. 

Changed  lower  absolute  latitude  boundary  of  low/mid-latitude  transition  region  from  30.  degrees  to 
34.  degrees. 

Changed  upper  absolute  latitude  boundary  of  pure  low-latitude  region  from  30.  degrees  to  34. 
degrees. 

Removed  comments  regarding  f0E  normalization  from  METHOD  comment  section  since  no  f0E 
normalization  is  done.  _ 


Removed  arguments  TCI  and  TC2  from  calls  to  routine  REGMOD  since  they  are  no  longer  used  by 
that  routine. 

Removed  local  variables  TCI  and  TC2  since  they  are  no  longer  used.  . _ 


Moved  the  call  to  routine  PHANTM  from  after  the  conversion  of  the  nominal  UT  from  hours  to 
seconds  to  before  the  call  to  routine  MIDLAT  so  that  phantom  data  does  not  influence  the  sunspot 
number  iteration. 

Removed  the  call  to  routine  INIT  that  accompanied  the  call  to  routine  PHANTM  since  it  is  no  longer 
needed. 

Updated  the  version  number  and  version  date. 


Changed  PARAMETER  MOPM1  from  1 1  to  14  for  new  LLF  parameterization. _ 


Removed  INCLUDE  statement  for  INCLUDE  file  tomid.inc  since  it  is  not  needed. _ 


Removed  arguments  EE,  EF,  FE,  FF  from  the  call  to  routine  DO_ADJ  since  they  are  no  longer 

needed  by  that  routine. _ _ _ _ 

Modified  adjustment  algorithm  to  use  a  unified  multiplicative  scaling  that  varies  smoothly  across 
altitude  boundaries. 

Removed  input  parameters  NMHME,  NMHMF2,  NOHME,  and  NOHMF2  since  they  are  no  longer 

used. _ _ _ _ _ 

Corrected  calculation  of  topside  half-width  by  changing  difference  HMH-ALTF(J)  to  ALTF(J)- 
IIMII. 
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MEMORANDUM 


DATE: 

29-March-2001 

TO: 

Files 

FROM: 

Rob  Daniell 

RE: 

Changes  to  PRISM  1. 7c  for  PRISM  1. 7d 

The  changes  to  PRISM  1.7c  for  PRISM  1.7d  focus  on  correcting  problems  in  both  the  high  latitude  and  midlatitude 

real-time  adjustment  algorithms,  bringing  the  use  of  TEC  data  in  line  with  PRISM  1.5  documentation,  and  providing 

a  logical  switch  to  control  whether  or  not  “phantom  stations”  are  used.  The  changes  are  summarized  as  follows: 

1.  BUG  FIX:  Dwight  Decker  discovered  that  when  PRISM  was  driven  with  TEC  data  alone,  no  phantom  stations 
were  assigned  even  though  there  were  large  regions  with  no  data.  This  problem  has  been  corrected  so  that 
phantom  stations  are  correctly  assigned  regardless  of  what  combination  of  ionosonde  and  TEC  data  are  ingested. 
The  changes  were  made  in  file  “phantom.for”.  (However,  note  item  2  below.) 

2.  A  compile  time  logical  parameter  PHANTOM  has  been  added  to  PRISM  to  control  whether  “phantom  stations” 
are  generated  to  force  die  model  to  relax  to  climatology  in  regions  where  data  is  absent.  For  operational  use, 
this  parameter  should  be  set  to  .FALSE,  to  suppress  the  generation  of  phantom  stations.  This  change  was  made 
in  the  file  “prism.for”. 

3.  BUG  FIX:  An  error  in  the  call  to  SUBROUTINE  OVAL  from  SUBROUTINE  FMODEL  was  corrected.  In 
addition  a  programming  errot  in  SUBROUTINE  OVAL  itself  was  also  corrected.  These  errors  were  identified 
by  Dwight  Decker  and  caused  problems  under  certain  conditions  in  the  auroral  oval  corrections.  The  changes 
were  made  in  file  “fmodel.for”. 

4.  BUG  FIX:  Several  lines  of  code  used  to  adjust  the  high  latitude  E-layer  was  “commented  out”  in  Version  1.7c 
without  any  accompanying  documentation.  This  caused  errors  in  the  high  latitude  adjustments  under  some 
circumstances.  The  lines  were  “uncommented”  in  Version  1.7d,  and  this  seems  to  have  corrected  these  errors. 
This  change  was  made  in  the  file  “init.for”. 

5.  The  time  limits  on  DMSP  data  were  changed  to  conform  to  PRISM  documentation.  PRISM  now  accepts  DMSP 
data  for  all  times  within  ±2  hours  of  the  nominal  time  of  the  run.  This  change  was  made  in  file  “rtdlimitinc”. 

6.  BUGFIX:  An  error  in  the  interpolation  of  EOF  coefficients  in  UT  was  corrected.  The  error  caused  incorrect 
calculation  of  the  interpolated  coefficients  when  data  from  more  than  one  time  are  ingested.  This  change  was 
made  in  the  file  “gen.for”. 

7.  The  coding  that  sets  a  floor  for  the  electron  density  was  corrected  to  ensure  that  the  floor  is  always  enforced. 
This  change  was  made  in  the  file  “gen.for”. 

8.  BUG  FIX:  The  threshold  electron  flux  and  proton  flux  for  identification  of  the  auroral  oval  boundaries  was 
corrected  to  conform  to  PRISM  documentation.  This  problem  was  discovered  by  Dwight  Decker.  This  change 
was  made  in  the  file  “getdat.for” 

9.  The  transition  altitude  between  the  E-layer  and  F-layer  in  the  EDP  adjustment  algorithm  has  been  lowered.  This 
should  prevent  the  formation  of  false  F-layer  peaks  between  the  true  F-peak  and  the  E-peak  when  the  E-layer 
adjustment  is  much  larger  than  the  F-layer  adjustment.  This  change  was  made  in  the  file  “rta.for”. 
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10.  The  algorithm  for  the  usage  of  TEC  data  has  been  modified  to  bring  it  into  conformance  with  PRISM  1.5 
documentation  (PL-TR-95-2061).  When  the  DISS  site  (with  valid  /F2  data)  nearest  to  a  particular  TEC  datum 
is  within  the  “decorrelation  distance”  DCL  (defined  as  a  PARAMETER  in  “prism.for”),  then  the  TEC  data  is 
used  to  calculate  a  pseudo-SSIES  record  unless  at  least  one  real  SSIES  datum  is  also  present  within  a  distance 
DCL  of  the  TEC  datum,  in  which  case  the  TEC  data  is  discarded.  Otherwise,  the  TEC  data  is  used  to  calculate  a 
pseudo-DISS  record.  Although  derived  from  the  SEND  version  of  PRISM,  the  implementation  differs 
somewhat.  The  EDP  parameter  corrections  fields  are  calculated  with  five  consecutive  passes  through  the  real¬ 
time  data.  (1)  The  first  pass  uses  DISS  data  to  calculate  F-  and  F-layer  height  corrections.  (2)  The  second  pass 

uses  DISS  data  to  calculate  E-  and  F-layer  peak  density  corrections.  (3)  The  third  pass  simply  locates  DMSP  r 

SSIES  data  for  later  comparisons  with  the  locations  of  DISS  and  IMS  data.  (4)  The  fourth  pass  uses  TEC  data 
to  calculate  pseudo-DISS  corrections  (if  the  nearest  DISS  site  is  beyond  the  “decorrelation  length”)  or  pseudo- 
SSIES  corrections  (if  the  nearest  DISS  site  is  within  the  “decorrelation  length”  and  there  is  no  SSIES  data  also 
within  the  “decorrelation  length”).  (5)  Finally,  the  fifth  pass  uses  the  DMSP  SSIES  data  to  calculate  topside 
corrections  for  the  EDP.  While  not  completely  self  consistent,  this  should  cause  the  correction  fields  to  more 
closely  reproduce  the  TEC  data  at  the  TEC  data  locations.  These  changes  were  made  in  the  files  “hlim.for”, 
“midlatfor”,  and  “newfitinc”  as  well  as  several  other  subroutines  that  call  SUBROUTINE  REGMOD. 

1 1.  The  parameter  limiting  the  amount  of  high  latitude  data  that  can  be  ingested  by  PRISM  was  increased  from  8000 
to  20000  by  Dwight  Decker.  This  change  was  made  in  the  file  “hlism.for”. 

12.  BUG  FIX:  When  By  vanishes  exactly,  PRISM  would  halt.  This  has  been  changed  so  that  By= 0  is  treated  as 
B}>0.  This  problem  was  identified  by  SMC  Det  11  personnel  during  testing  of  PRISM  with  actual  data.  This 
change  was  made  in  the  file  “iojutil.for”. 

13.  A  compile  time  logical  parameter  DIAGNOS  has  been  added  to  PRISM  1.7d  to  control  certain  diagnostic  output 
that  is  useful  for  debugging  and  code  development.  While  this  diagnostic  output  has  been  designed  not  to 
interfere  with  operational  use  of  PRISM,  it  does  alter  the  standard  output.  Therefore,  we  strongly  recommend 
that  for  normal  operations,  DIAGNOS  should  be  set  to  .FALSE,  in  the  file  “prism.for”.  When  DIAGNOS  is 
TRUE,  there  are  two  separate  kinds  of  diagnostic  output  that  are  produced: 

a.  The  region  code  (CAP,  OVAL,  TROUGH,  MIDLAT)  is  included  in  the  standard  output.  In  gridded  output, 
it  is  added  to  the  end  of  the  line  containing  the  geographic  and  geomagnetic  coordinates  of  each  grid  point. 

In  station  output,  its  location  depends  on  the  user’s  choice  of  output.  If  the  user  has  selected  EDP  data  only, 
the  region  code  appears  immediately  after  the  magnetic  local  time  (as  for  gridded  output).  If  the  user  has 
selected  profile  parameter  data  only  or  both  EDP  and  profile  parameters,  the  region  code  appears 
immediately  after  the  profile  parameters. 

b.  PRISM  writes  a  file  named  “boundary.txt”  containing  the  parameters  that  determine  the  boundaries  between 
the  trough,  the  amoral  oval,  and  the  polar  cap.  An  example  file  is  shown  below.  The  printed  parameters 
correspond  to  the  parameters  described  in  PL-TR-95-2061  (PRISM:  A  PARAMETERIZED  REAL-TIME 
IONOSPHERIC  SPECIFICATION  MODEL,  VERSION  1.5)  on  pages  25-26  and  in  Appendix  D.  The  first 
line  of  the  file  includes  the  name  of  the  gridded  output  file  if  it  exits.  (If  there  is  no  gridded  output,  then  the 
first  line  simple  reads  “Boundary  Parameters”.  The  second  line  is  the  actual  Kp  as  provided  in  the  user 
input.  The  third  line  gives  the  effective  Kp  for  the  amoral  F-layer.  The  remainder  of  the  file  gives  the 


Boundary 

Parameters 

for  "prism_21.out"  1 

3.0000 

-  actual  Kp 

3.0000 

«  "effective"  Kp  for 

auroral  E- layer 

3.0000 

3.0000 

3.0000  = 

"effective"  Kp  for  F-layer  boundaries 

0.0000 

0.0000 

0.0000  = 

THETA0 (0:2) 

-10.5000 

0.0000 

0.0000  = 

AB (0:2) 

11.5000 

0.0000 

0.0000  = 

BB { 0 : 2 ) 

3.8800 

0.0000 

3.0000  - 

CB (0:2) 

2.7300 

0.0000 

0.0000  = 

PB { 0 : 2 ) 

-10.5000 

2.7000 

0.8000  = 

CAB (0:2) 

o.nooo 

0.2670 

-0.2670  = 

VAB (0:2) 

24.4000 

20.9000 

13.4000  = 

CTH (0:2) 

2.1200 

1.7000 

1.7000  = 

VTH (0:2) 

11.5000 

2.6333 

=  CBB (0:1) 

0.0000 

-0.0833 

=  VBB (0:1) 

0.0000 

0.0051 

=  QBB (0:1) 

* 
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parameters  that  define  the  F-layer  boundaries  as  used  in  SUBROUTINE  FMODEL.  They  are  described  in 
the  header  to  the  file  “regionjb.inc”  which  defines  the  boundary  parameters  contained  in  die  block  common 
BNDS.  This  header  is  displayed  below. 


C  Declarations  for  high  latitude  region  boundaries 

C  As  of  6  April  1990  the  trough  equatorward  boundary  is  a  distorted  circle 
C  of  the  form 
C 

C  THETA ( PHI )  =  THETAO ( Kp )  +  AB*EXP ( - (ABS ( (PHI -BB) /CB) apb) ) 

C 

C  where  THETA  is  magnetic  co- latitude  (degrees),  PHI  is  magnetic  local  time 
C  (hours)  .  THETAO  definitely  depends  on  Kp;  AB,  BB,  CB,  and  PB  may  or  may 
C  not.  The  equatorward  auroral  oval  boundary  is  an  undistorted  circle  whose 
C  center  is  displaced  from  the  magnetic  pole.  The  polar  cap  boundary  is  also 
C  a  displaced  circle.  For  these  THETAO  is  the  radius  (degrees)  ,  AB  is  the  CGM 
C  co- latitude  of  the  center,  BB  is  the  MLT  of  the  center. 

C 

C  At  low  magnetic  activity,  the  precipitation  oval  and  the  trough  overlap 
C  during  midmoming  (0800-1000)  magnetic  local  times.  This  may  (or  may  not) 

C  be  an  indication  of  the  absence  of  a  morning  trough  at  low  activity 
C  levels.  CB  (1)  is  a  minimum  displacement  of  the  oval  boundary  from  the 
C  trough  boundary,  i.e.,  a  minimum  width  of  the  trough.  Likewise,  CB{2)  is  a 
C  minimum  width  of  the  oval.  PB(1)  and  PB(2)  are  currently  unused. 

C 

C  CB  (0)  =  the  local  time  scale  size  of  the  distortion: 

C 

C  CB(1:2)  are  the  minimum  trough  and  oval  widths 

C  PB(0)  is  the  exponent: 

C  PB(1:2)  are  currently  unused 

C 

C  THETAO ( 0 )  =  24.4  +  2.12*KP(0) 

C  THETAO ( 1 )  =  20.9  +  1.7*KP(1) 

C  THETAO (2)  “  13.4  +  1.7*KP(2) 

C 

C  And  fix  the  amplitude  and  local  time  of  the  distortion  of  the  trough 

C  boundary: 

C  AB ( 0)  =  -10.5 

C  AB ( 1)  =  2.7  +  0 . 267*KP (1) 

C  AB (2 )  =  0.8  -  0 . 267*KP (2) 

C 

c 

REAL  THETAO  (0:2)  ,  AB  (0  :  2 )  , BB  (0 : 2)  , CB  (0 :  2 )  ,  PB  (0 : 2)  ,  CAB  (0 : 2 )  , 

1  VAB  (0:2)  ,  CTH  (0:2)  ,VTH(0:2)  ,CBB(0:1)  ,VBB(0:1)  ,QBB(0:1) 

C 

COMMON/BNDS /THETAO ,  AB ,  BB ,  CB ,  PB ,  CAB ,  VAB ,  CTH ,  VTH ,  CBB ,  VBB ,  QBB 


14.  A  number  of  mainly  cosmetic  changes  were  made  to  improve  readability  and  internal  documentation  of  the 
code,  to  provide  diagnostic  output  to  the  log  file,  to  improve  portability,  and  to  simplify  test  runs  during  code 
development.  These  changes  are  transparent  to  the  user. 

15.  BUG  FIX:  A  logic  error  in  the  DMSP  data  ingestion  routine  INDMSP  (in  file  “getdat.for”)  was  corrected.  If 
there  was  a  blank  line  or  an  erroneous  section  label  line,  the  routine  would  loop  through  the  data  file 
indefinitely.  Now,  if  a  line  cannot  be  identified  as  belonging  to  one  of  the  three  data  types  (drift,  in  situ  plasma, 
or  precipitating  particles)  the  line  is  skipped,  and  the  next  line  is  read.  Note  that  this  bug  would  have  no  effect 
on  properly  formatted  DMSP  data  files. 

16.  The  parameter  MDISS,  defined  in  the  file  “rtdlimitinc  ”  was  changed  to  MDISS  =  41000+MPIONO  to 
accommodate  the  potentially  large  number  of  pseudo-DISS  sites  to  be  generated  from  SSUSI  disk  images. 


The  table  on  the  next  page  describes  the  changes  that  I  made  to  PRISM  1.7c  to  produce  PRISM  1.7d. 
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Module 


Program  Unit 


phantom.for 


ftnodel.for 


fmodel.for 


getdat.for 


output. for 


output. for 


output. for 


output. for 


output. for 


output. for 


midlat.for 


hlism.for 


init.for 


io i  util. for 


rtdlimit.inc 


newfit.inc 


SUBROUTINE  PHANTOM 


SUBROUTINE  FMODEL 


SUBROUTINE  OVAL 


SUBROUTINE  INDMSP 


SUBROUTINE  OUTPUT 


SUBROUTINE  GRID_OUTPUT 


SUBROUTINE 

STATIONOUTPUT 


SUBROUTINE  WRITEJDATA 


SUBROUTINE  W_ST_DATA 


SUBROUTINE  WRTBND 


SUBROUTINE  DO_ADJ 


SUBROUTINE  GEN  E 


SUBROUTINE  GTCFF 


SUBROUTINE  REGMOD 


SUBROUTINE  MIDLAT 


SUBROUTINE  HLISM 


SUBROUTINE  INIT 


SUBROUTINE  UDETID 


N/A 


Description  of  Changes  _ 


(1)  Modified  the  way  SUBROUTINE  PHANTM  write  phantom  ionosonde  data  to  the  scratch  files 
that  PRISM  uses  to  organize  the  real-time  data  it  has  ingested.  This  corrected  a  problem  the  failure 
of  PRISM  to  store  phantom  station  data  when  IMS  data  was  ingested  in  the  absence  of  DISS  data. 

(2)  Added  the  logical  parameter  PHANTOM  to  control  whether  phantom  station  data  is  generated 
(PHANTOM  =  .TRUE.)  or  whether  phantom  stations  are  suppressed  (PHANTOM  =  .FALSE.).  The 
value  of  PHANTOM  is  printed  at  the  beginning  of  the  log  file. 

(3)  Added  the  logical  parameter  DIAGNOS  to  controls  whether  or  not  RCODE  is  appended  to  the 
standard  output  and  whether  or  not  SUBROUTINE  WRTBND  is  called  to  produce  the  file 
“boundary.txt”  containing  the  boundary  parameters  for  the  trough,  auroral  oval,  and  polar  cap 
boundaries.  The  value  of  DIAGNOS  is  printed  in  the  log  file  immediately  below  PHANTOM. 

(4)  Added  the  logical  parameter  LTECTOP  to  control  whether  or  not  IMS  data  that  is  “near”  DISS 

data  is  used  to  generate  pseudo-SSIES  data  or  pseudo-DISS  data.  The  real  parameter  DCL  specifies 
the  “TEC  decorrelation  length”  that  is  used  to  determine  whether  IMS  data  is  “near”  DISS  data.  The 
values  of  these  parameters  are  now  printed  at  the  beginning  of  the  log  file  immediately  after 
PHANTOM  and  DIAGNOS.  _ 


Modified  the  way  phantom  stations  are  assigned  so  that  it  is  done  correction  regardless  of  the 
combination  of  DISS  and  IMS  data.  _ 


Corrected  out  of  order  variables  in  the  call  to  OVAL 


Corrected  error  in  linear  formula.  (Note:  cubic  formula  should  be  reinstated  and  tested  when  time 
permits.)  _ 


Changed  PARAMETER  CRTEF  from  2.5  to  0.25  to  conform  to  documentation. 

Changed  PARAMETER  CRTPF  from  1.0  to  0.1  to  conform  to  documentation. 

Corrected  logic  error  in  SUBROUTINE  INDMSP  which  could  potentially  cause  the  routine  to  enter 
an  indefinite  loop  condition  in  the  case  of  certain  errors  in  the  format  of  the  DMSP  data  files. _ 


Added  DIAGNOS  to  the  argument  list  and  to  the  calls  to  GRID_OUTPUT  and  STATION__OUTPUT. 
Changed  STATUS  parameter  from  ‘NEW’  to  ‘UNKNOWN’  in  OPEN  statements  for  output  files. 

This  allows  PRISM  to  overwrite  old  output  files  rather  than  crashing.  _ 


(1)  Added  latitude  region  code  (CHARACTER*8  RCODE)  to  standard  output. 

(2)  Added  DIAGNOS  to  argument  list  to  detemine  whether  or  not  RCODE  is  appended  to  standard 
output. 

(3)  Added  TECTOP  to  call  to  REGMOD. 

(4)  Added  diagnostic  print  statements  to  permit  the  observation  of  the  profile  adjustment  process  at  a 

specific  grid  point.  (See  also  SUBROUTINE  DO_ADJ  in  “rta.fof”.)  These  print  statements  are  set 
off  bv  lines  of  “#####”  symbols.  _ 


Added  TECTOP  to  call  to  REGMOD.  Added  DIAGNOS  to  argument  list  and  to  call  to 
W_ST_DATA.  _ 


Added  RCODE  and  DIAGNOS  to  argument  list.  Added  logic  to  determine  whether  or  not  to  append 
RCODE  to  the  write  statement  that  writes  GLAT,  GLON,  etc,  depending  on  the  value  of  DIAGNOS. 


Added  RCODE  and  DIAGNOS  to  argument  list.  Added  logic  to  determine  whether  or  not  to  append 
RCODE  to  write  statements.  The  location  of  RCODE  in  the  station  output  depends  on  the  type  of 
output  selected  by  the  user.  _ 


Created  this  subroutine  (called  from  PROGRAM  PRISM)  to  write  latitude  boundary  parameters  to 
file  “boundary.txt”.  _ 


(1)  Modified  the  way  the  transition  height  between  the  E-  and  F-layers  is  calculated  to  prevent  (or 
reduce  the  chances  of)  the  formation  of  false  FI -layers  below  the  true  F2-layer  peak. 

(2)  Added  diagnostic  write  statements  (controlled  by  a  logical  variable  set  in  GRID_OUTPUT  and  ■ 

passed  in  common  block  RTAFL)  to  allow  the  profile  adjustment  process  to  be  monitored. _ 


Modified  final  test  on  DEN  to  insure  that  densities  are  no  smaller  than  DENMIN _ 


Corrected  error  in  logic  that  caused  incorrect  interpolation  in  UT  when  near  a  midnight  crossing  and 
data  from  more  than  one  UT  has  been  ingested.  _ 


Modified  argument  list  and  code  to  allow  for  a  calculation  of  the  “topside  TEC”  (the  TEC  above  the 
peak  density)  needed  for  the  corrected  TEC  ingestion  algorithm. 


Modified  code  to  allow  for  the  corrected  TEC  ingestion  algorithm. 


Increased  MMODEL  from  8000  to  20000,  _ _ _ 


Restored  call  to  LDITER  for  E-layer  height  adjustment  (lines  296-299 


Corrected  handling  of  B 


Modified  limits  on  DMSP  data  so  that  PRISM  accepts  DEMSP  data  for  all  times  within  ±2  hours  of 
the  nominal  time  of  the  run.  Modified  the  limits  on  DISS  data  to  accommodate  the  potentially  large 
number  of  pseudo-DISS  sites  to  be  generated  from  SSUSI  image  data.  _ 


Changed  the  parameters  MNTD  and  MHTD  to  equal  MISP+MIMS  to  reflect  the  possible  addition  of 
pseudo-SSIES  records  from  TEC  data  collocated  with  DISS  data. 
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M  E-MO  R  A  N  D*  UM 


DATE: 

9-January-2002 

TO: 

Files 

FROM: 

Rob  Daniell 

RE: 

Changes  to  PRISM  1. 7dfor  PRISM  1. 7e 

The  changes  to  PRISM  1.7d  for  PRISM  1.7e  are  principally  designed  to  correct  a  problem  with  the  assimilation  of 
TEC  data  in  the  presence  of  SSIES  data.  However,  other  changes  were  intended  to  bring  PRISM  1.7  into 
compliance  with  PRISM  1.5  documentation,  particularly  with  regard  to  the  assimilation  of  TEC  data,  and  the 
remaining  changes  were  designed  to  make  the  coding  easier  to  follow. 

1.  The  data  ingestion  algorithm  in  SUBROUTINE  MIDLAT  was  modified.  It  still  makes  four  passes  through  the 
data,  and  the  first  two  passes  continue  to  search  for  h„F2  data  and  foF2  data,  respectively.  Pass  3  now  searches 
for  TEC  data  that  is  farther  than  the  “decorrelation  length”  (DCL,  a  PARAMETER  in  the  main  program)  from 
the  nearest  DISS  site  and  uses  it  to  calculate  a  pseudo-DISS  data  record.  Pass  4  now  searches  for  both  SSIES 
and  TEC  data.  TEC  data  that  is  farther  than  DCL  from  the  nearest  SSIES  datum  is  used  to  calculate  a  pseudo- 
SSJES  record.  TEC  data  that  lies  closer  than  DCL  to  an  SSIES  datum  is  ignored.  In  connection  with  this 
change,  SUBROUTINE  PRO_COR  and  SUBROUTINE  COR_PRO  were  completely  rewritten  to  conform  to 
Version  1.5  documentation  regarding  the  topside  profile  correction  algorithm.  These  changes  were  made  in  the 
file  “midlat.for”  and  “rta.for”. 

2.  The  variable  TOPHT  has  been  added  to  the  calls  to  SUBROUTINE  PARAM  and  SUBROUTINE  FMODEL  in 
SUBROUTINE  REGMOD  and  SUBROUTINE  RTA.  Previously,  the  variable  HTOP  did  double  duty  as 
TOPHT  (the  altitude  of  the  topside  density  correction)  on  input  and  HTOP  (the  topside  scale  height)  on  output. 
Now  the  variables  are  separate  entities,  resulting  in  cleaner,  easier  to  understand  code.  These  changes  were 
made  in  the  files  “hlim.for”,  “fmodel-for”,  “param.for”,  and  “rta.for”. 

3.  BUG  FIX:  Dwight  Decker  modified  the  loop  over  DMSP  data  types  in  SUBROUTINE  INDMSP  to  correct  a 
problem  reading  data  from  multiple  satellites.  He  also  added  some  information  messages  to  the  log  file  from 
this  subroutine.  The  changes  were  made  in  file  “getdat.for”. 


The  table  on  the  next  page  describes  the  changes  made  to  PRISM  1.7d  to  produce  PRISM  1.7e. 


» 
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Module 

Program  Unit 

Description  of  Changes 

midlat.for 

SUBROUTINE  MIDLAT 

Modified  the  data  ingestion  algorithm.  The  revised  algorithm  makes  four  passes  through  the  DMSP 
data.  The  first  pass  searches  for  AmF2  data  from  DISS  sites  and  calculates  height  corrections.  The 
second  pass  searches  for  fJFi  data  from  DISS  sites  and  calculates  peak  density  corrections.  The  third 
pass  searches  for  TEC  data  that  are  farther  than  the  “decorrelation  length”  (PARAMETER  DCL 
specified  in  PROGRAM  PRISM)  from  the  nearest /*F2  datum  and  calculates  a  peak  density 
correction.  The  fourth  pass  searches  for  both  SSIES  and  TEC  data.  Topside  density  and  scale  height 
(where  applicable)  corrections  are  calculated  from  the  SSIES  data  and  from  TEC  data  that  are  farther 
than  DCL.  TEC  data  that  are  closer  than  DCL  to  at  least  one  SSIES  datum  are  ignored. 

midlat.for 

SUBROUTINE  PRO_COR 

Rewrote  subroutine  to  conform  to  Version  1 .5  documentation  regarding  topside  corrections. 

hlim.for 

SUBROUTINE  REGMOD 

Added  TOPHT  to  argument  list  so  that  HTOP  would  not  change  meaning  between  input  and  output 

fmodel.for 

SUBROUTINE  FMODEL 

Added  ZTOP  to  argument  list  and  changed  HTSAV  to  HTOP 

param.for 

SUBROUTINE  PARAM 

Added  ZTOP  to  argument  list  so  that  HTOP  would  not  change  meaning  between  input  and  output 

rta.for 

SUBROUTINE  RTA 

(1)  Added  ZTOP  to  argument  list  so  that  HTOP  would  not  change  meaning  between  input  and 
output.  Added  ZTOP  to  call  to  SUBROUTINE  CORJPRO. 

(2)  Changed  topside  correction  algorithm  so  that  it  is  invoked  only  if  DNTP  or  DHTP  is  nonzero. 

rta.for 

SUBROUTINE  COR_PRO 

Rewrote  subroutine  to  conform  to  Version  1.5  documentation 

getdat.for 

SUBROUTINE  INDMSP 

Modified  loop  over  data  types  (SSIES  ION  DRIFT,  SSIES  IN  SITU  PLASMA,  and  SSJ/4  DATA) 
according  to  recommendation  of  Dwight  Decker  to  correct  problems  reading  data  from  multiple 
satellites. 

> 
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